Cold polar molecules in 2D traps: 
Tailoring interactions with external fields for novel quantum phases 
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We discuss techniques to engineer effective long-range interactions between polar molecules using 
external static electric and microwave fields. We consider a setup where molecules are trapped in a 
two-dimensional pancake geometry by a far-off-resonance optical trap, which ensures the stability 
of the dipolar collisions. We detail how to modify the shape and the strength of the long-range part 
of interaction potentials, which can be utilized to realize interesting quantum phases in the context 
of cold molecular gases. 



I. INTRODUCTION 



The realization of Bosc Einstein condensates and quan- 
tum degenerate Fermi gases with cold atoms has been 
one of the highlights of experimental atomic physics dur- 
ing the last decade and in view of recent progress in 
preparing cold molecules we expect a similarly spectac- 
ular development for molecular ensembles in the co ming 
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|2(X |21| . The outstanding features of the physics of cold 
atomic and molecular gases are the microscopic knowl- 
edge of the many-body Hamiltonians, as realized in the 
experiments, combined with the possibility to control and 
tune system parameters via external fields. Examples are 
the trapping of atoms and molecules with magnetic, elec- 
tric and optical traps, allowing for the formation of quan- 
tum gases in ID, 2D and 3D geometries, and the tuning of 
contact inter-particle interactions by varying the scatter- 
ing length via Feshbach resonances [IE HE . This control 
is the key for the experimental realization of fundamen- 
tal quantum phases, as illustrated by the superfluid-Mott 
insulator quantum phase transition with bosonic atoms 
in an optical lattice 1241 ] , an d the BEC-BCS crossover in 
atomic Fermi gases [gl [20, HE HE [29J 

As discussed in our recent work [3f3 |. polar molecules 
prepared in the electronic and vibrational ground state 
offer new possibilities to control inter-particle interac- 
tions. In fact, effective interactions with a given potential 
shape can be engineered under conditions of tight 2D 
confinement, by applying static (DC) and microwave 
(AC) fields. The engineered potentials can display both 
repulsive and/or attractive character. This control of 
the interactions - in combination with low-dimensional 
trapping - opens the way to realizing novel quantum 
phases and quantum phase transitions. As an example 
Ref. [30| discusses a quantum phase transition from a 
superfluid to a self-assembled crystal for a gas of polar 
molecules in the strongly interacting limit, where the 
stability of the collision processes is guaranteed by the 



confinement in a 2D geometry It is the purpose of the 
present paper to present in some detail the molecular 
aspects behind this engineering of effective two-body 
interactions. 
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FIG. 1: (color online) System setup: Polar molecules are 
trapped in the (x,y)-plane by an optical lattice made of two 
counter-propagating laser beams with wavevectors ±kL = 
iki^e z (arrows on the top/bottom). The dipoles dj are 
aligned in the z-direction by a DC electric field Edc = Edc&z 
(arrow on the left). An AC (microwave) field Eac is indi- 
cated (arrow on the right). Inset: Definition of polar (■&) and 
azimuthal (ip) angles for the relative orientation of the inter- 
molecular collision axis r with respect to a space-fixed frame 
with axes {e x ,e y ,e z }. 

The interaction potential between atoms, in particular 
Alkali atoms in their electronic ground state, is domi- 
nated at large distances by an attractive Cq /r 6 potential. 
In the many body Hamiltonian for a dilute quantum gas 
this gives rise to an effective two-body short range inter- 
action in the form of a contact interaction with a scat- 
tering length a s . Polar molecules have strong permanent 
electric dipole moments in their electronic-vibrational 
ground state manifold, and pairs of molecules aligned by 
external DC or AC electric fields will interact via (com- 
paratively strong) dipolc-dipolc interactions with char- 
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acteristic long-range 1/r 3 dependence 0, 
These dipole-dipole interactions will be attractive or re- 
pulsive, depending on the relative orientation of the 
dipoles. 

The alignment of the dipoles corresponds to the dress- 
ing of the lowest energy excitations of the internal molec- 
ular degrees of freedom, which are related to rotations of 
the molecule. The rotational dynamics can be manip- 
ulated using external electric DC and AC (microwave) 
fields. This dressing of rotational states by external fields 
together with the dipole-dipole interaction forms the ba- 
sis to shape the effective molecular interactions. 




FIG. 2: (color online) Qualitative sketch of effective 2D poten- 
tials Vejp(p) for polar molecules confined in a 2D (pancake) 
geometry. Here, p = r sin i?(cos ip, sin ip) is the 2D coordinate 
in the plane z = and p = rsini?, (see inset of Fig. [TJ. Solid 
line: Repulsive dipolar potential Vjp(p) = C^/p 3 induced by 
a DC electric field. Dash-dotted line: "Step-like" potential in- 
duced by a single AC (microwave) field and a weak DC field. 
Dashed line: Attractive potential induced by the combination 
of several AC (microwave) fields and a weak DC field. Here, 
the potentials V^g (p) and the separation p are given in arbi- 
trary units. For the "step-like" case (dash-dotted line) p = 1 
corresponds to the Condon point rc of Sec. IIII CI 

One example discussed in Ref. [3(| deals with polar 
molecules confined in a 2D (pancake) trap (see Fig. []}. 
The molecular dipoles are aligned perpendicular to the 
plane by a DC field. Thus, the effective 2D interac- 
tions are repulsive and long range l/JjP(p) = C^/p 3 , with 
p = r sin?9(cos</?,sin<^) the 2D coordinate in the plane 
2 = and p = rsin$ (see Fig. [2J solid line). The in- 
teraction strength C3 is proportional to the square of the 
induced dipole moment for the dressed rotational ground 
state. Depending on the interaction strength, we find 
the appearance of a crystalline phase, and an associated 
quantum melting to a superfluid phase as a function of 
the square of the induced dipole- moment. The corre- 
sponding phase diagram is discussed in Ref. [3(|, and it 
is reproduced in Fig. [3] (see also Ref. [H|). 

In the present work we present in detail the micro- 
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FIG. 3: (color online) Sketch of the phase diagram for a ho- 
mogeneous 2D system of polar molecules interacting via the 
effective 2D repulsive potential V 2 g{p) = C-j/p 3 . T is the 
temperature in units of Td = Cs/ksa 3 , with a the aver- 
age inter-particle distance and &b the Boltzmann constant. 
The symbol ra = E lnt /E^\ n = Czm/h 2 a is the interaction 
(-Bint = Cs/a 3 ) to kinetic energy (-Ekm = h 2 /ma 2 ) ratio. A 
crystalline phase appears for large ratios > t-qm and small 
temperatures T < T m . The critical ratio tqm ~ 18 ±4 for the 
quantum melting to a superfluid phase has been determined in 
Ref. [3(|, while the classical melting temperature T m (dashed 
line) to a normal gas phase has been calculated in Ref. [36l [. 
The finite-temperature superfluid to normal fluid phase tran- 
sition is of the Berezinskii-Kosterlitz-Thouless type [37l . l38l ] 
and it appears below the upper bound T KT = ttH 2 /2k B ma 
(dashed-dotted line) . The crossover to an unstable regime for 
small repulsion and finite confinement in the z-direction (see 
Fig. [1} is indicated by a hatched region (see text, Sec. III). 

scopic molecular theory underlying this engineering of 
the interaction potential for trapped polar molecules of- 
fered by DC and AC microwave fields. We focus both 
on potentials which arc repulsive 1/r 3 (DC field) and 
on potentials which have a marked "step-like" character, 
that is, the character of the repulsive potentials varies 
considerably in a small region of space (an AC plus a 
DC field). Three example cases of effective 2D potentials 
Vtfp(p) are shown in Fig. (5) The use of multi-chromatic 
AC fields can lead to the realization of interesting po- 
tentials (for example the attractive potential of Fig. rj|), 
however, in this work we focus on monochromatic AC 
fields only. 

In all cases, the derivation of the effective 2D inter- 
actions proceeds in two steps: First, we derive a set 
of Born-Oppenhcimcr (BO) potentials by diagonalizing 
the Hamiltonian for the relative motion of two particles 
for fixed molecular positions. Within an adiabatic 
approximation, the corresponding eigenvalues play the 
role of an effective 3D interaction potential. Second, 
we obtain an effective 2D dynamics by integrating-out 
the fast transverse motion of the molecules along the 
direction of the tight parabolic confinement. 
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The paper is organized as follows: In Sec. |TT] we dis- 
cuss the Hamiltonian for a single rotating polar molecule 
dressed by DC and AC (microwave) fields under con- 
ditions of strong optical confinement. The collisions of 
two polar molecules are considered in Sec. IIIII After re- 
viewing the molecular collisions in the absence of exter- 
nal fields, in Sec. IIII Al we consider the case of interac- 
tions in the presence of a DC electric field. In particular, 
the stabilizing effects of a parabolic potential confining 
the particles to a 2D plane are analyzed in Sects. HTl A 31 
and IIII A 41 while the effective 2D interaction potential 
y c 2 ff D (p) = C 3 /p 3 is derived in Sec. IIII A 51 The in- 
teractions in the presence of an AC field are studied 
in Sec. IIII Bl In the absence of external confinement, 
this case is analogous to the 3D optical shielding de- 
veloped in the context of ultracold collisions of neutral 
atoms H(| As in the latter, we find a strong 

dependence of the 3D shielding potential on the polar- 
ization of the AC field. The presence of "holes" in the 
3D shielding potential for certain polarizations and of 
several degeneracies in the two-particle spectrum for all 
polarizations render the pure AC-ficld case less appealing 
for realizing stable collisional setups in 2D. In fact, both 
the former and the latter processes open loss channels 
for the ground-state interaction. In Sec. IIII O we analyze 
the interactions in combined DC and AC fields, and we 
show that the DC field helps to greatly suppress the pres- 
ence of possible loss channels at large distances, while an 
additional harmonic confinement along z avoids populat- 
ing the regions of space where "holes" analogous to those 
of Sec. IIII Bl occur. Thus, by introducing a tight optical 
confinement in the z-direction, in this case it is possible 
to realize stable two-dimensional collisional setups. Two- 
dimensional inter-particle interactions can be designed, 
whose character varies markedly between long and short 
distances, allowing for much greater flexibility in tuning 
by external fields than the pure DC case of Sec. IIII Al 



II. MOLECULAR HAMILTONIANS 

The purpose of this section, which forms the basis 
of discussion in the following sections, is to review the 
single-molecule rotational spectroscopy. In particular 
we are interested in the rotational excitations of cold 
v 2S+1 A(v) spin-less (S = 0) polar molecules in their elec- 
tronic (y = 0) and vibrational (v = 0) ground-state, with 
zero-projection (A = 0) of the total angular momentum 
on the internuclear axis [42|, |43[ . The spectroscopic nota- 
tion for the electronic-vibrational ground-state of these 
molecules is X 1 S(0). Moreover, we are interested in ma- 
nipulating the rotational states of these molecules using 
DC and AC electric fields and in confining the particles 
using a (optical) far-off-resonance trap (FORT). The ap- 
plication of these external fields will serve as a key ele- 
ment to engineer effective interaction potentials between 
the molecules. 

Our goal in this section is to derive a low energy effec- 



tive Hamiltonian for the external motion and internal ro- 
tational excitations of a single molecule in its electronic- 
vibrational ground state of the form 

H(t) = ^ + H Iot + H uc + H AC (t) + tf pt(r). 
Ira 

In the last equation, p 2 /2m is the kinetic energy for the 
center-of-mass motion of a molecule of mass m, while 
H lot accounts for the rotational degrees of freedom. The 
terms i^DCj HAc(t) and H opt (r) refer to the interaction 
with electric DC and AC (microwave) fields and to the 
optical trapping of the molecule in the ground electronic- 
vibrational manifold, respectively. 



A. Rotational excitations of E molecules 

We consider spin-less polar molecules with E electronic 
ground-states in their electronic-vibrational ground 
state, X 1 E(u = 0). The low-energy internal excita- 
tions correspond to the rotation of the internuclear axis 
of the molecules with total internal angular momentum 
J [42|, HH, [44|. The corresponding Hamiltonian H rot is 
the one of a rigid spherical rotor [42j 

H Iot =BJ 2 . (1) 

Here B is the rotational constant for the electronic- 
vibrational ground state, which is of the order of 
B ~ h 10 GHz J45[. We denote the energy eigenstates 
of Eq. (fTJ) by \J,M), where J is then quantum number 
associated with the total internal angular momentum 
and M is the quantum number associated with its 
projection onto a space-fixed quantization axis. The 
excitation spectrum is Ej — BJ(J + 1), which is 
anharmonic. Each J- level is (2 J + l)-fold degenerate. 

A polar molecule has an electric dipole moment, d, 
which couples its internal rotational levels. This dipole 
moment gives rise to the dipolc-dipolc interaction be- 
tween two molecules. For E-molecules the dipole opera- 
tor is along the internuclear axis e a b, i.e. d = de a b. Here, 
d is the "permanent" dipole moment of a molecule in its 
electronic-vibrational ground-state. 

The spherical components of the dipole operator on 
a space-fixed spherical basis {e_i,eo,ei}, with e q= o = 
e 2 and e±i = T( e x ± i e y)/v / 2, are given by d q = e q ■ 
d = dC$\e,<t>) ,whcre cf } (<?,</>) arc the unnormalizcd 
spherical harmonics and 9 (</>) is the polar (azimuthal) 
angle for the orientation of the molecule in the space- 
fixed frame [H, 0, [44| , respectively. We note that for 
a spherically-symmetric system, e.g. in the absence of 
external fields, the eigenstates of the rotor have no net 
dipolc-momcnt, (J, M\d\J,M) = 0. On the other hand, 
the component d q couples the rotational states | J, M) 
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and | J ± 1 , M + q) according to 
(J ±1,M + q\d q \J,M) = d(J,M; l,q\J±l,M + q) x 
x(J,0; 1,0|J± 1,0), 



2J+ 1 



2(J± 1) + 1' 

where (Ji, Mi; J 2 , M 2 \J, M) arc the Clebsch-Gordan- 
cocfficicnts. 

In the following we are interested in the interaction of 
the molecules with an external DC electric field along 
e 2 , Edc = Euc^o, and with AC microwave fields with 
cither linear polarization (q = 0) or circular polarization 
(q = ±1) relative to e z , Eac(£) = EA_ce~ lut e q + cc. 
These fields couple to a molecule via the electric dipole 
interaction, 

Hoc = -d-E DC = -d Euc, (2a) 
H AC (t) = -d-E AC (i) = -d q E acb-^ +h.c, (2b) 

which try to align the molecule along the field, while 
competing with its rotation, as [J 2 ,<i g ] ^ 0. 



B. Coupling of rotational states by DC and AC 
electric fields. 
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FIG. 4: (color online) Solid lines: Energies Ej : m (left) and 
states |0j,m) (right) of Eq. (j4| with J = 0, 1, for a molecule in 
a weak DC electric field Edc = Euc^o with f3 = dEuc/B <C 
1. The DC-field-induced splitting H5 and the average energy 
separation to are M = 3B/3 2 /20 and to = 2B + Bf3 2 /6, 
respectively. Dashed and dotted lines: Energy levels for a 
molecule in combined DC and AC fields (The AC-Stark shifts 
of the dressed states are not shown). Dashed line: The AC 
field is monochromatic, with frequency u>, linear polarization 
q — 0, and detuning A = lj — (Q + 25/3) > 0. Dotted lines: 
Schematics of energy levels for an AC-field with polarization 
q = ±1 and frequency u' ^ uj. 



1. Coupling to a DC electric field 

The effects of a DC electric field, Edc, on a single 
polar molecule are: (a) To split the (2J + l)-fold 
degeneracy in the rotor spectrum, and (b) to align the 
molecule along the direction of the field, which amounts 
to inducing a finite dipole moment in each rotational 
state. 

We choose the direction of the DC-field as the quanti- 
zation axis, Edc = -Ebc e o- Then, the internal Hamilto- 
nian is that of a rigid spherical pendulum [Z^, [46| 



H — -ff r ot + Hoc — BJ — doE 



DC, 



(3) 



which conserves the projection of the angular momentum 
J on the quantization axis, i.e. M is a good quantum 
number. The energy eigenvalues and eigenstates of 
Eq. are labeled as Ej^u and |0j,m), respectively. 

We are interested in weak fields, Edc <C B/d, where 
the effects of the electric field are a quadratic DC Stark 
shift of the rotational energy levels and a finite induced 
dipole moment along the axis of the field in each ro- 
tational state. For a typical rotational constant, B ~ 
h 10 GHz, and a dipole-moment d ~ 9 Debye this cor- 
responds to considering DC fields (much) weaker than 
B/d ~ 2 kV/cm. To lowest order in = dE DC /B the 
energy eigenvalues and eigenstates are [H, |4o| 



Ej,m/B 



1 ; 2 (2J-l)(2J + 3) V ' 



2 yjW+l) 



p vV + i) 2 -m 2 

"2^+1)3(2.7+1) 



\J + 1,M). (4b) 



Thus, the ground state energy is shifted downwards by 
i?o,o = —Bf3 2 /6, while the energies of the lowest excited 
states are split by 



hS = E lt0 - E h±1 = 3B/3 2 /20, 



(5) 





(00,0 do 0O,o) 


(d/3/3)(l - 7/3 2 /360) 


gi 


(01, ±1 \do 0i, ±i) 


(d/3/10)/(l-3/3 2 /5600) 


92 


(0i,o|do|0i,o) 


-(d/3/5)(l - 19/3 2 /350) 


fo 


(0i,o do I0o, o) 


{d/VS)(l -43/3 2 /360) 


h 


(0i, ±i \d±i 0o, o> 


(d/v / 3)(l- 49/3 2 /1440) 


h 


(0i,o Mti 01,±l) 


(3d/J/20)(l + ll/3 2 /1400) 



TABLE I: Permanent (g n ) and transition (/„) dipole mo- 
ments of the 4 states belonging to the rotational (J = 0, 1)- 
manifolds, in the presence of a weak polarizing DC field, 
Edc = -Edc^o- Here, f3 = dEuc/B <§C 1 is the ratio of the 
electrostatic energy and rotational constant, and the dipole 
moments are given up to order third order in j3. 
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see solid lines in Fig. |4j The average energy separation 
of the (J = 0) and (J = l)-manifolds is 

1 

m = (Ei,m- E .o)/3 = 2B + B(3 2 /6. (6) 

M=-l 

The induced dipole moments to lowest order in (3 are 

3M 2 /J(J+ 1) - 1 



>j,m u <f>J,M 



dp- 



(2J-l)(2J + 3) 



■e . 



This equation shows that the ground state acquires a fi- 
nite dipole moment go = (0o,oMo|0o,o) — dp/ 3 along the 
field axis, while the lowest excited states acquire a dipole 
moment (<^i,Af |do|^i,M> = d(3M 2 - l)/3/10. For later 
convenience, perturbative values in the small parameter 
P of the transition and induced dipole moments are re- 
ported in Tablc|T]for the four single-particle states \^>j,m) 
with \M\ < J < 1. The transition and induced dipole 
moments are labeled as /„ and g n , respectively. 



2. Coupling to an AC electric field 

Similar to the case of a DC electric field, the basic 
effect of an AC electric field Eac on a single molecule 
is to polarize it by dressing its energy levels. The 
characteristic time-dependence of the AC field allows 
for: (a) Addressing individual rotational transitions by 
applying one or several non-interfering microwave fields 
(multi-mode field); (b) Realizing dressing fields that 
can be not only linearly, but also circularly polarized. 
In this work, we consider the case of a single AC 
microwave field with polarization q and frequency uj, 
EacC 1 ", t) = E A c(r)e~ lult e q + ex., and derive the dressed 
energy-levels for a molecule in the field. For the sake of 
generality - and for later convenience - we consider the 
case where the AC field is superimposed to a weak DC 
field, which provides for a splitting of the first excited 
( J = 1)— manifold, as shown above. 

Given a polarization q, the frequency uj is chosen close 
to the transition from the ground-state to one state of the 
(J = l)-manifold, |<^>o,o) l^i,?), i- c - w ~ ^ + 5(2/3 — 
q 2 ), where the states |</>j,m) are those of Eq. The 
corresponding wavelength is of the order of centimeters, 
which largely exceeds the size of our system and therefore 
one can neglect the position dependence of the microwave 
field, i.e. recoil effects, E A c( r ) ~ -^ac- The electric 
dipole interaction of Eq. (|2b|) reads 



H AG (t) 



-d„EAce~ 



h.c. 



(7) 



The Rabi frequency and 
are f2 = EAc(<t>i,q\d q \<f>ofl)/K 
A = u - (E hq - E 0>0 )/h = uj 
respectively (see Fig. 0]). 



the detuning A 
= E Ac f\ q \/h and 
p + 6(2/3 - q 2 )}, 



In Sec. IIIII we consider a specific setup where the 
AC field has linear polarization, q = 0. Here we il- 
lustrate how to obtain the dressed energy levels of a 
molecule in this field by diagonalizing the Hamiltonian 
H = H Toi +Hnc+H\c(t) in a Floquet picture. First, we 
expand the Hamiltonian on the basis \4>j : m), which diag- 
onalizes the time- independent part of H as iJ ro t +-Hdc = 
J2j m \4>J,m)Ej,m(4>J,m\- Then, we consider the effect of 
the AC field driving the (|0o,o) ^ |</>i,o))-transition with 
Rabi-frequency £1 = /o-Eac ~ dE ac / \/3h and detun- 
ing A = uj - (E h0 - E 0>Q )/h = 25(1 + 2p 2 /15)/h. 

A transformation to the Floquet picture is obtained by 
expanding the time-dependent wave- function in a Fourier 
series in the AC frequency to. After applying a rotating 
wave approximation, i.e. keeping only the energy con- 
serving terms, we obtain the time-independent Hamilto- 
nian H, which describes the coupled two-level system in 
the basis {|0o,o), \<Pi,o)} as 



H 



Q 
CI A 



Eq,o- 



The corresponding dressed energy eigenvalues of H for 
the ground state and excited state (minus one photon 
energy huj) are given by 



-E'0,0 — -Eo,o — 



HA 
hA 



hA 
~Y 
hA 



An 2 

'A 2 



hn 2 



i 



4^ 2 

'A 2 



+ - 



-HA- 



W 2 



respectively. We note that the AC field induces an AC- 
Stark shift ss ±fi,il 2 /A, on the ground and the excited 
state, respectively. Thus, the shift depends on the de- 
tuning A, and in particular on its sign, and on the Rabi- 
frequency n. 



C. Optical trap 

An essential ingredient of our setup is the tight 
confinement of the molecules in a 2D-plane. This is 
realized for example by a far-off-resonant optical trap 
(see Fig. []}. The latter drive far-off- resonance transi- 
tions from XS(0) to the electronically excited states, 
i/A(v). The goal of this section is to obtain the resulting 
trapping potentials for the lowest rotational excitations, 
J = 0,1. 

A detailed discussion of the complex nature of molec- 
ular electronic excitations [42], |43[ is beyond the scope of 
the present discussion. For a detailed treatment of an 

We consider here a 



example case we refer to Ref. [47 [. 
simple model, where the fine and hyperfine interactions 
are neglected. Then the basic molecular structure is ob- 
tained as follows: In the adiabatic approximation one 
diagonalizes the Hamiltonian for the electrons and the 
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two nuclei as a function of the internuclear separation 
r a f,, thus obtaining a set of Born-Oppcnhcimcr (BO) po- 
tentials. E u ^(r a b). Here v is the main electronic quantum 
number, while A denotes the quantum number associated 
with the operator for the total angular momentum com- 
ponent of the molecules along the internuclear axis, e a f J. 
The latter gives rise to a large splitting of the electronic 
manifolds ~ A„|A| 2 , where A v = h 2 / I e is the inverse of 
the (small) moment of inertia of the electrons [42j . Then, 
the vibration of the nuclei in the BO-potentials yields a 
series of bound-states v = 0, 1, 2, . . . with energy E v \r v \ . 

Deep optical traps are obtained as follows: The laser, 
E opt (r, t) = E opt (r)e~ lWL * + c.c, drives the electronic 
transitions to the lowest excited states, labeled A and 
B, with frequency wl tuned near the minima of the BO- 
potcntials. Since spontaneous emission in the excited 
states is typically a few MHz, deep traps on the order 
of one MHz with a negligible inelastic scattering rate (~ 
a few Hz) require detunings on the order of hundreds of 
GHz from the vibrational resonances. Since these detun- 
ing are much larger than B, one can neglect the rotational 
structure in the electronic ground and excited states in 
deriving the optical potential. The effective interaction 
of the molecules with the off-resonant laser-field is thus 
described by 

H opt (r) = E opt (r)*-d( WL )-E opt (r), (8) 

with the dynamic polarizability-tensor 

a(ui h ) = a\\(u L )e' <g> e + a±(u L ) ^ (-l) A e^ <8> e'_ A 

A=±l 



a_i_ 



(o> L ) ^2 (-1)% <g> e_ 9 + [a||(w L ) - aj_(u L )] 



g=-l 



(9) 

Here {e'_ 1; eg, e' +1 } denotes a body-fixed spherical basis 
with e = e ab = S g (-l) 9 C_g (0, <t>)e q being the inter- 
nuclear axis, and a||(wL) [aj_(wL)] is the dynamic polar- 
izability at frequency ojl in the direction parallel (per- 
pendicular) to the internuclear axis [13, EH ■ The parallel 
and perpendicular components are given by the S — £ 
(AA = 0) and S - n (AA = ±1) transitions, respec- 
tively, and read 



«nK) = EE 



|di/S(n)-XE(0)| 



, (10a) 



) = EE 



\d v n(v)-xz(o) | 2 



(10b) 



Here d,/A(ii)-xs(o) denotes the transition dipole moment 
from the XE(0) to uA(v), and the sum over ± accounts 
for the near-resonant and typically far off- resonant terms. 
From Eq. (|10p we see that the anisotropy in the dynamic 



polarizabilitics, o;j_(wl) — a ||( w L), is due both to the dif- 
ferent dipole-moments and to the large splitting of the 
excited vYi{v) and vli{v) states. 

In our setup we consider a pair of circularly po- 
larized counter-propagating laser beams, E opt (r) = 
E ov t cos(£;L£)e_|_, with wave-vectors ±1cl = iuj-^eo/c 
along z, trapping the molecules in the x — y plane (see 
Fig. [T|). From Eq. ([5]) we obtain the following Hamilto- 
nian for the optical trapping [48[ 

i?opt(r) = a a (uj h )\E opt \ 2 cos 2 (k L z)C { °\e^)+ 

+a 2 (co L )\E opt \ 2 cos 2 (k h z)C^ 2) (6, 0), (11) 

where a (cj L ) = [ot\\(wi,) + 2aj_(uj L )]/3 and a 2 (^L) = 
[a||(wL) — a^(a>L)]/3. The first term in Eq. (|lip . pro- 
portional to c^ ] {e,4>) = 1, gives an overall shift, which 
is common to all the rotational states. The second term 
is responsible for tensor shifts, which split the excited 
rotational states according to |M|, as 



(j,M\ci 2 \e,4>)\j,M) 



J(J + 1) - 3A/ 2 
(2J- l)(2J + 3)' 



(12) 



Typical depths of optical lattices are of the order of 
< A 1 MHz, and thus much smaller B. Therefore we 
may neglect the far-off resonant Raman coupling between 
different J manifolds, i.e. J <-> J ± 2. 

We consider tight optical traps, such that the 
molecule in the ground state are strongly confined 
at one potential minimum of (<^o,o|-ffopt(r)|^o,o) = 
ao{uJL)\E op t\ 2 cos 2 (kLz). For a light field which is 
(far) red detuned from the electronic excited states, i.e. 
hu>L -C E U M V ) ~ Exs(o), the dynamic polarizabilitics 
a II (o>l) and o;_l(wl) are negative and the trapping poten- 
tial for the ground-state is attractive, since cco^l) < 0- 
We assume the molecule to be strongly confined near the 
field anti-node z = 0. Then the optical trapping is essen- 
tially given by a tight harmonic trap 

ffopt(r) « \a (u L )\\E opt \ 2 (-l + k 2 z 2 ) + 

+a 2 (io h )\E opt \ 2 (l - k 2 z 2 )C^ 2 \d, <p).(13) 

From the last expression we see that the tensor-shifts 
induce a position-dependent splitting for the excited ro- 
tational manifolds, which at z = is analogous to that 
induced by a DC field, but it has a strong modulation 
in space. The tensor-shifts are thus seen as position and 
state-dependent potentials, and the last term in H opt (r) 
is (in principle) unwanted for our purposes, since it gives 
rise to different trapping frequencies uj± for the ground 
and excited states. 

However, we note that by applying a second laser, 
t (r, i), of frequency cu' L with wavevector and po- 
larization e^, one can eliminate the state-dependent 
potentials - up to a position independent splitting of 
the excited states. Given the large number and vari- 
ety of available excited electronic-vibrational states sev- 
eral choices are possible. One choice is, e.g., to ap- 
ply an additional laser with the same polarization ej^ 
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as the first laser, i.e. = e+, but having a node at 
3 = and being blue detuned from the electronic transi- 
tions, i.e. E' opt (r,t) — E' t sm(k' h z)e~ llAj ^ t e + + c.c. with 
< » (#„A(0) - E xm )/H for vk = AH, BY. This 
induces an additional state dependent optical trapping 
potential given by H' opt {r) = \E' opt \ sin 2 (A£*)[a (wL) + 
a 2 (uj' L )CQ '(9,<p)]. Tuning the laser-frequency cu' L with 
respect to the vibrational resonances one can force both 
c*_l(<^l) an d a ||( w L) to be positive, see Eq. (JTUJ). The ad- 
ditional trapping potentials are zero at the node z = 0, in 
particular for the ground-state the trapping potential is 
repulsive (thus enhancing the trapping given by the first 
laser), while the excited-state position-dependent trap- 
ping oc z 2 c^ 2) (e,(f>) of Eq. CCD can be compensated for 

by tuning the strength of the second laser, E' t . 

The parabolic trapping potential for our setup is then 
given by 

ff op t(r) = ^mulz 2 - V + V 2 C^(9, 0), (14) 

where the first term is a state-independent har- 
monic trapping along e z at frequency u>± = 
[2\a (u; L )\\E opt k L \ 2 / m _ + 2a Q {uj' l )\E' opt k< l ] 2 / m] 1 1 2 , 
the second terms gives an overall Stark-shift, 
Vo = |ao( w L)||£'opt| 2 j and the last term is a split- 
ting of the excited rotational states, J > 0, which is 
independent of the position z, V2 = a^^iJI-E'optl 2 - 

Concluding, the Hamiltonian for a single molecule is 

H(t) = |- + l -mulz 2 -V + V 2 C^ (0, 0)+ 

+BJ 2 - d a E DC - (d q E Ac e~ wt + h.c.) . (15) 



Cq ip) = e q ■ e r are unnormalized spherical hamon- 
ics with $ and tp polar and azimuthal angles relating the 
orientation of e r with respect to a space-fixed frame e q , 
respectively. The terms d qv j = e q ■ dj are the spherical 
components of the projection of the dipole operator of 
molecule j onto the space- fixed frame e q . 

In the absence of external fields -Edc = ^ac = 0, 
the interaction of the two molecules in their rotational 
ground state is determined by the van-der-Waals attrac- 
tion Kdw ~ C6:o/?* 6 with Cq : o ~ —d 4 /6B. This expres- 
sion for the interaction potential is valid outside of the 
molecular core region r > tb = (d 2 /B) 1 / 3 , where tb 
defines the characteristic length where the dipole-dipole 
interaction becomes comparable to the splittings of the 
rotational levels, see below. In the following we show that 
it is possible to induce and design interaction potentials 
which arc long-range, by dressing the interactions with 
appropriately chosen static and/or microwave fields. In 
fact, the combination of the latter with low-dimensional 
trapping allows to engineer effective potentials whose 
strength and shape can be both tuned. The derivation 
of the effective interactions proceeds in two steps: (i) We 
derive a set of Born-Oppcnheimcr (BO) potentials by 
first separating Eq. (fl"6)) into center-of-mass and relative 
coordinates, and diagonalizing the Hamiltonian for the 
relative motion for fixed molecular positions. Within an 
adiabatic approximation, the corresponding eigenvalues 
play the role of an effective 3D interaction potential in a 
given state manifold dressed by the external field, (ii) We 
eliminate the motional degrees of freedom in the tightly 
confined z direction to obtain an effective 2D dynamics 
with interaction V 2 ^(p). In the following we consider 
the cases of a static field and a microwave field, coupling 
the lowest rotor states. 



III. TWO MOLECULES 

We consider the interactions of two polar molecules 
j = 1 , 2 confined to the x — y plane by a tight harmonic 
trapping potential of frequency directed along z. The 
interaction of the two molecules at a distance r = r 2 — 
T\ = re r is described by the Hamiltonian 



(16) 



where Hj(t) is the single- molecule Hamiltonian Eq. (fT"5)) . 
and Vdd(r) is the dipole-dipole interaction 



^dd(r) 



di ■ d 2 - 3 (di • e r ) (e r • d 2 



(17) 



Here, d^- is the dipole operator of the molecule j, and 
e r -dj is its projection onto the collision axis e r . The pro- 
jection reads e r ■ dj = Y^t=-i(~ l) 9 C_g($i f)d q -j, where 



Effective interactions in the presence of a DC 
electric field 



In this section we consider the collisions of two ground- 
state molecules in the presence of a DC electric field, 
Edc = ^Dceo- The Hamiltonian Eq. p^|) now reads 



H = 



2 

E 

2 

E 

3=1 L 



^ + \mJiz) + BJ 2 - Eocdw 
P? 1 2 2 



V dd {r) 



tfint(r), 



(18) 



where do-.j = eo • dj and iJ; n t(r) is the internal Hamil- 
tonian including the dipole-dipole interaction, iJj nt (r) = 
J2j [-^Jj — E~Dcdo;j] + Vdd(r), respectively. In this sec- 
tion we are interested in ground-state collisions, and thus 
for convenience we set V% = in Eq. (|T5]) . that is, we ne- 
glect possible tensor-shifts in the excited-state energies 
of each molecule. We can further rewrite Eq. (fl~8"|) by 
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splitting H in center of mass and relative coordinates as 
H = H com + H Iel with 

p2 

H com = —+mujiZ 2 , (19a) 
4m 

ffrel = — + ]mcjlz 2 + H int (r). (19b) 

TO 4 

Here, R = (ri + r2)/2 and P = pi + P2 are the center 
of mass coordinate and momentum of the two molecules, 
while r = r2 — ri and p = (p2 — Pi)/2 are the relative 
coordinate and momentum, respectively. Equations (|19|) 
show that the dipole-dipole interaction couples the inter- 
nal degrees of freedom to the relative motion, while the 
latter and the harmonic motion of the center of mass re- 
main decoupled. Thus, the non-trivial system's dynamics 
is entirely determined by H TC \. In the following we focus 
our discussion on this term. 

As explained above, in the spirit of the BO- 
approximation we can obtain effective interaction poten- 
tials for the collision of the two particles by diagonal- 
izing f/rei for fixed particle positions and zero kinetic 
energy. In the adiabatic approximation, the resulting 
eigenvalues are energy-surfaces which act as effective po- 
tentials in each state manifold. Here we first analyze the 
case of collisions in the absence of external fields, that 
is E DC = w_l = 0, Sec. MI A 11 Then, in Sec. MI A 21 we 
add a static electric field of small strength -Edc 
The effects of finite trapping uj± ^ are treated in the 
following section Sec. MI A 31 for the most relevant case 
of ground-state collisions. The stability of ground-state 
collisions is investigated in Sec. MI A~4| The effective two- 
dimensional potential for ground state collision is derived 
in Sec. MI A 51 

1. Collisions in the absence of external fields 

In the absence of external fields (-Edc = <*>j_ = 0) and 
for zero kinetic energy, diagonalizing H IC \ amounts to di- 
agonalizing if; nt (r) as a function of r, 

2 

H int (r) = ]T BJ 2 + V dd (r) = £ |$„(r))£„(r)($„(r)|, 

j—1 n 

(20) 

where E n (r) and |$ n (r)) are the n th -adiabatic energy 
eigenvalues and two-particle eigenfunctions, respectively, 
and n is a collective index for a set of quantum num- 
bers to be specified below. Each eigenvalue E n (r) 
plays the role of an effective interaction in a given 
state manifold dressed by the external field. At in- 
finite separations of the molecules, the eigenfunctions 
|$! 0) (r)) = |$i 0) (tf,^)) = |$ n (r ->oo,0,¥>)> are sym- 
metrized products of the (rotated) single-particle eigen- 
states \J j} Mj)j = e ~ ivJ ^e~ MJ y^ \ Jj, Mj)j, which are 
independent of the distance r. For finite r the two- 
particle cigenstatcs are superposition of several single- 
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FIG. 5: (color online) BO-potentials _Bj ; y ;CT (r) as a function 
of the distance r, for two molecules interacting in the absence 
of external fields, Euc = u± = /3 = (see text and Tab. Hi). 
Here, J = Ji + J2, Y and a are the total number of rota- 
tional excitations shared by the two molecules, the quantum 
number associated with the projection of the total internal 
angular momentum onto the collision axis, = e r • (Ji + J2), 
and the permutation symmetry under the exchange of the par- 
ticles, respectively. The solid and dashed curves correspond 
to symmetric (a — +) and antisymmetric (a — — ) eigen- 
states, respectively. Each potential energy surface is labeled 
by the corresponding energy and eigenstate (see Tab. 111)1 . 
Here, rs = (d 2 / B) 1 ^ 3 , with d the permanent dipole moment 
and B the rotational constant of each molecule, respectively. 
Note that the n- and the A-states are doubly degenerate. 



particle states, which are mixed by the dipole-dipole in- 
teraction Vdd(r). 

A few eigenvalues E n (r) of Eq. ([20]) are plotted as 
a function of r in Fig. \5\ Figure O shows that the 
energy spectrum behaves quite differently for r < tb 
and r > tb, tb = (cP/B) 1 / 3 . In fact, for r < rs a large 
number of level crossings and anticrossings occurs, which 
make the fulfillment of the adiabatic approximation 
generally impossible. The region r < rs is the molecular 
core region. In the following we focus on the region 
r > rs, where the lowest-energy eigenvalues group into 
well defined manifolds, which are approximately spaced 
by an energy ps 2B. For ground-state collisions, the 
adiabatic approximation is here trivially fulfilled. 
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Since we are interested in ground-state collisions, we 
restrict our discussion to the Jj = and Jj = 1 manifolds 
of each molecule, which amounts to take into account 16 
rotational two-particle states. The corresponding eigen- 
values E n (r) are clearly distinguishable in Fig. [5] in the 
region r > rg. The manifolds are approximately split 
by 2B, according to the number of rotational excitations 
Ji + J2 shared by the two molecules. The two-particle 
energy eigenstates and eigenpotentials can be classified 
according to the following symmetries of -ffint(r): (a) 
The projection of the total internal angular momentum 
along the collision axis, e r ■ (Ji + J2), is conserved and 
associated with a quantum number Y; (b) The Hamil- 
tonian is invariant under the exchange of the two parti- 
cles, which is associated with the permutation symmetry 
a = ± under the exchange of the two particles. This 
implying that symmetric (antisymmetric) states couple 
to symmetric (antisymmetric) states only; And (c) the 
parity p = a(—l) Jl+J ' 2 is conserved. The spectroscopic 
notations labeling the eigenstates, |$ n (r)), and poten- 
tials, E n (r), with n = (J = Ji + J2\Y\o) in Fig. [5] are 
explained in the caption of Tab. [TTJ 

Analytic results for the energy eigenvalues and eigen- 



n 


J 


Y 


a 


Y, 


E (0) 


C 3 ;„/d 2 


C 6; „ x 6B/d 4 











+ 










-1 


1 


1 





+ 


s u 


2B 


-2/3 


-22/45 


2,3 


1 


±1 




n g 


2B 


-1/3 


-19/45 


4,5 


1 


±1 


+ 


n u 


2B 


+1/3 


-19/45 


6 


1 







s g 


2B 


+2/3 


-22/45 


7,8 


2 


o± 


+ 


s g 


4B 





-(48=f39V3)/50 


9,10 


2 


±1 




n u 


4B 





-39/20 


11, 12 


2 


±2 


+ 


A g 


4B 





-24/25 


13, 14 


2 


±1 


+ 


n g 


4B 





-51/25 


15 


2 









4B 





-6/25 



TABLE II: Perturbative expansions of the effective poten- 
tials, E n = Ej.y,a, for the three lowest-energy manifolds, 
J = J1+J2 = 0,1,2 with Jj — 0,1, for interactions in 
the absence of external fields Euc = Eac = — 0. First 
column: index n = 0, 1,2, . . . labeling the energy potentials 
and states. Second column: total number of rotational ex- 
citations shared by the two molecules J = Ji + Ji. Third 
column: the quantum number Y associated with the projec- 
tion of the total internal angular momentum along the col- 
lision axis, e r • (Ji + J2), which is denoted by E, II, A for 
Y = 0,1,2, respectively. For n = 7, 8 the subindex ± indi- 
cates the presence of two states, which are split by the Van- 
der-Waals interaction. Fourth column: permutation symme- 
try a = ±. Fifth column: spectroscopic notation with parity 
p — <r(— 1) j denoted by "g" (gerade) for p — +1 and "u" 
(ungerade) for p — —1, respectively. Sixth column: Asymp- 
totic energy E^ = E n (r — > 00). Seventh column: dipole- 
dipole coefficient C3 ;n . Last column: Van-der-Waals coeffi- 
cient Ce- n . Perturbative energy eigenvalues are expressed in 
the form: E n (r) = E n a) + C 3; „/r 3 + C 6; „/r 6 . 



states of -£/i nt (r) for large enough inter-particle dis- 
tances r can be derived using a perturbative expan- 
sion in Vdd(r)/i?. Our results for the energy eigen- 
values E n (r) = E Jl+ j 2;Y -A r ) of B T IJ j2 j + Vdd(r) are 
summarized in Tab. [TTJ There, the asymptotic energy 
En^ = E n (r — > 00), the dipolc-dipole coefficient C^ n 
and the Van-der-Waals coefficient C§. n are reported, so 
that the perturbative expression for E n (r) takes the form 

E n (r) = E ( n } + C 3; „/r 3 + C 6; „/r 6 . The Table shows 
that the ground-state energy £h( r ) = -Eb;0;+( r ) is shifted 
downwards by an amount i?o,o,+( r ) = — d 4 /6Br e , which 
is the usual Van der Waals shift due to off-resonant 
dipolc-dipole interactions. The first excited manifold, 
(Ji + J2 = 1), consists of 6 states, of which 3 are sym- 
metric and 3 are antisymmetric. These states are split by 
the resonant dipolc-dipole interaction according to their 
angular momentum along the collision axis, \Y\ = 0, 1 
and a = ±, as reported in Table ITT1 Finally, the second 
excited manifold, (Ji + J2 = 2), consists of 9 states, of 
which 6 are symmetric, Y = 0±, ±1, ±2 with a = +, and 
3 arc antisymmetric, Y = 0, ±1 with a = — . 



2. Collisions in a DC field: Effective 3-D interaction 

We now turn to study the collision of the two molecules 
in the presence of a weak static electric field applied in the 
^-direction but in the absence of optical trapping, that is, 
E = i?DC e o with EVjc <C B/d and oj± = 0. As explained 
in Sec. Ill B 11 the effects of a DC electric field on each 
molecule are to partially split the (2 J+l)-fold degeneracy 
in the rotor spectrum (the modulus of the projection M is 
conserved) , and to align the molecule along the direction 
of the field, which amounts to inducing a finite dipolc 
moment ji^j^ ,Mj Mo;j \4>Jj,M ■) 3 m eacn rotational state. 

Analogous to the discussion above, the effec- 
tive interaction potentials for the collision of the 
two particles can be obtained in the adiabatic 
approximation by diagonalizing the Hamiltonian 
H Icl = H int {r) = E 3 2 =i [Btf - Encd . tj ] + F dd (r) = 
J2 n |$„(r))-E„(r)(<I>„(r)|, where now the asymptotic 
energy eigenstates |<&i ^(t9, if)) are symmetrized products 
of the single-particle states Ifij^M^j of Eq. (|4]). The 
quantity n = (J; M; a) is the collective quantum number 
labeling the eigenvalues E n (r), with J = Ji + J%, 
M = \Mi \ + |M 2 |, and a = ±. We note that, because 
of the presence of the DC field, here J is a simple label 
for the various energy manifolds, and not a quantum 
number. The energies of the eigenvalues E n (r) and the 
associated eigenvectors are tabulated in Tab. IIII1 

Similar to the zero-field discussion, in the weak field 
limit (3 <C 1 and for r > tb we expect the eigenval- 
ues of i/i n t(r) to group into manifolds, which are ap- 
proximately separated by the rotational spacing IB. On 
the other hand, because of the finite induced dipolc mo- 
ments j (4>j u m \do-j\4>j ,m )j, for the two molecules can 
now interact resonantly via the dipole-dipole interaction 
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0.5 1 1.5 r/r B 1 10 r/r B 1 r* 10 r/r B 



FIG. 6: (color online) BO-potentials E.j-,M- a (r, •&) for two molecules colliding in the presence of a DC field, with j3 = dEuc/B — 
1/5 and {J\M\a) the quantum numbers of Tab. IIIII The solid and dashed curves correspond to symmetric (a = +) and 
antisymmetric (o = — ) eigenstates, respectively, (a): BO-potentials for the 16 lowest-energy eigenstates En(r,-&). The 
molecular-core region is identified as the region r < rs = (d 2 / B) 1 ^ 3 , while for r ^> tb the eigestates group into manifolds 
separated by one quantum of rotational excitation 2B. (b) and (e): Blow-ups of the first-excited energy manifold of panel (a) in 
the region r > tb for ■& = ir/2 and $ = 0, respectively. Note the electric-field-induced splitting HS = 3B f3 2 /20 (see Sec. II B 1). 
The distance rs where the dipole-dipole interaction becomes comparable to hS is rs = (d 2 /hS) 1 ^ 3 . (c) and (e): Blow-ups of 
the ground-state potential Eofi;+(r, $) of panel (a) in the region r > tb for •& — n/2 and $ = 0, respectively. The distance r*, 
cf. Eq. (|22|l . where the dipole-dipole interaction becomes comparable to the Van-der-Waals attraction is indicated. Note the 
repulsive (attractive) character of the potential for ■& — n/2 (i9 = 0) and r > r*. 



Vdd(r) in each state manifold. This has important con- 
sequences for ground-state collisions. In fact, the new 
effective ground-state potential Eq(t) derived in pertur- 
bation theory in Vdd (r) / B reads 

^3:0 2 o\ ^6:0 



V™(r)^E Q (r) 



(l -3 cos 2 1?) 



(21) 



where a constant term 2i?o,o = —f3 2 B/3 due single- 
particle DC Stark-shifts has been neglected. The con- 
stants Cs-o ~ d 2 (3 2 /9 and C@ ; o ~ —d 4 /6B are the dipo- 
lar and Van der Waals coefficients for the ground-state, 
respectively (see Tab. IIIII) . Equation (|21|) is valid for 
t 3> re, and it shows that for distances r»r 4 with 



r* = 



2IC, 



C: 



:i:0 



1/3 



3d 2 



1/3 



(22) 



the dipole-dipole interaction dominates over the Van 
der Waals attractive potential, and V^P(r) ~ C3 ; o(l — 
3cos 2 i?)/r- 3 [3(J. In fact, the potential has a local maxi- 
mum in the plane z = r cos $ = at the position r = r*, 
where the dipole-dipole and Van der Waals interactions 
become comparable. The height of this maximum is 



Q 



3;0 



4C. 



6:0 | 



B(3 A 
~54~ 



(23) 



and the curvature along z is d 2 V(r 



-6C 3; o/^ 
quency 



0) 



-muj 2 /2, which defines a characteristic fre- 



12Q 



3:0 



1/2 



(24) 



to be used below. The latter has a strong dependence 
/3 8/3 = (dE DC /B) 8/3 on the applied electric field. 

We notice that if it were possible to confine the colli- 
sional dynamics to the (z = 0)-planc, purely repulsive 
long-range interactions with a characteristic dipolar 
spatial dependence ~ 1/r 3 could be attained. In the 
following sections, we analyze the conditions for realizing 
sufficiently strong confinements to the [z = 0)-planc 
by employing a tight harmonic optical trap in the 
z-direction. 

Figure [5] shows the eigenvalues E n (r) as a function 
of the interparticle distance r, for j3 = 1/5. The vec- 
tor r is expressed in spherical coordinates r = (r, i?, <f), 
with z = r cos^. Figure [6{ a) shows the different be- 
havior of the energy spectrum for r < r B and r > r Bl 
analogous to the zero-field case. Even at finite j3 we 
can clearly distinguish the molecular core region r < r B 
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E±e T ^l0i,o;0i,±i)/2 + perm. 


14 


2 


1+ 




2h(ZJ + 8/6) 




-3(27- 19T)/100 


E± e = F 1 * , |0i,o;0i,±i)/2 + perm. 
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2 





+ 


2h(uj + 2S/3) 


5 2 2 T 


-3(34-14T-T 2 )/100 


0i, o; 0i, o) 



TABLE III: Perturbative expressions for the 16 lowest-energy BO-potentials E n (r) = £i' J) + C 3 ; n h n (-ff)/r 3 + C 6; n(tf)/r 6 of two 
molecules interacting in the presence of a DC electric field Edc = (B/3/d)eo. First column: The collective quantum number 
n = (J — J\ + J 2 ; AI = A/i + |71/2|;<t — ±), labeling the eigenstates E n (r). Second column: The number J = J\ + J 2 of 
rotational excitations shared by the two molecules. Because of the presence of the DC field, parity is not conserved and J is a 
simple index that labels the various energy manifolds for r ^> rB,rs- Third column: The quantum number M = \M\ \ + \M 2 \. 
The additional subindex fj, for M > labels superposition of states with the same (Y;M;a), which depend on the azimuthal 
angle ip (see eigth column). Fourth column: The quantum number a = + (cr = — ) denoting symmetric (antisymmetric) states 
under permutation of the two molecules. Fifth column: Asymptotic energies for infinite separation. The quantities hS and 
ftZJ are defined in Eq. ((5} and Eq. JBJl, respectively. Sixth column: The C 3 - n coefficient and the angular dependence The 
dipole moments g n and /„ are defined in Table [T] while the angular distribution T is T = 1 — 3cos 2 i?. Seventh column: The 
contributions to the C 6; n(tf) coefficient up to order 0({3 2 ). The values A x , A 2 and A 3 are Ai = 40(2-T-T 2 )(/ /i+/ 2 go) 2 /d 4 /3 2 , 
A 2 = 33 + 6T-T 2 /2, A 3 = 13(1 + T)/ cos f, respectively. Here, £ is defined by the relation tan£ = (14 — T)(2 + T) /26(1 + T) . 
Last column: Eigenstates (1?, tp)) = |$ n (r — > 00, ■d.ip)} valid at infinite separation. Here, "perm." denotes the permuted 
state, e.g. |0i,2;03,4) — * |03,4;0i,2). 



where the adiabatic approximation breaks down. In this 
plot, the continuous and dashed lines correspond to the 
cases $ = 7r/2 and # = 0, respectively, which are almost 
indistinguishable on the scale of the graph. Figures[BJb,c) 
and Figs. [6jd,e) are blow-ups of the two lowest-energy 
manifolds of Fig. [SJa), for -d = ir/2 and i3 = 0, respec- 
tively. Different from the zero-field case, Fig. [|3b) and 
Fig. E^d) show that the excited state manifold with one 
quantum of rotation (J\ + J2 = 1) is asymptotically split 
into two sub-manifolds. This separation corresponds to 
the electric-field-induced splitting of the Jj = 1 manifold 
of each molecule, and it is thus given by HS = ZB0 2 / '20 
of Eq. ([5]). More importantly, Fig. [6](c) and Fig. [6]Je) 
show that the effective ground-state potential has a very 
different character for the cases d — ir/2 and d = 0, 
respectively. In fact, for i9 = tt/2, corresponding to colli- 
sions in the (z = 0)-plane [see Fig.^c)], the potential is 
repulsive and decaying at large distances as 1/r 3 in agree- 
ment with the discussion above. On the other hand, for 
t> = [see Fig. EJe)] the potential is purely attractive, 
with dipolar character. As mentioned above, in the next 
section we show that the probability to sample this at- 



tractive part of the potential during the collision can be 
largely suppressed in the case ui± 7^ 0, for a sufficiently 
tight transverse trapping. 

3. Parabolic confinement 

The presence of a finite trapping potential of frequency 
uj± in the z-direction provides for a position-dependent 
energy shift of Eq. (f2"Tj) . The new potential reads 

V(r) ee F e | D (r) + -mulz 2 

- ^(l-3cos^) + ^ + ^iz 2 .(25) 

As noted before, for z = the repulsive dipole-dipole 
interaction dominates over the attractive Van der Waals 
at distances r 3> given in Eq. (f2"2")l . In addition, for 
luj_ > the harmonic potential confines the particle's 
motion in the z direction. Thus, the combination of 
the dipole-dipole interaction and of the harmonic con- 
finement yields a repulsive potential which provides for 
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a three-dimensional barrier separating the long-distance 
from the short-distance regime. If the collisional energy 
is much smaller than this barrier, the particle's motion is 
confined to the long-distance region, where the potential 
is purely repulsive. 
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FIG. 7: (color online) Contour plot of the effective poten- 
tial V(p, z) of Eq. (|25|l . for two polar molecules interacting 
in the presence of a DC field (3 > 0, and a confining har- 
monic potential in the ^-direction, with trapping frequency 
lu ± = w c /10, where lu c = (12C 3; o/m^) 1/2 of Eq. (gtj and 
rv = (2|C6;o|/C 3 ;o) 1/3 of Eq. (HU. The contour lines are 
shown for V{p,z)/V+ > 0, with V* = Bf3 4 /54. Darker regions 
represent stronger repulsive interactions. The combination of 
the dipole-dipole interactions induced by the DC field and of 
the harmonic confinement leads to realizing a 3D repulsive po- 
tential. The repulsion due to the dipole-dipole interaction and 
of the harmonic confinement is distinguishable at z ~ and 
zjti, ~ ±7, respectively. Two saddle points (circles) located 
at (p±,±z±) separate the long-distance region where the po- 
tential is repulsive ~ 1/r 3 from the attractive short-distance 
region. The gradients of the potential are indicated by dash- 
dotted lines. The thick dashed line indicates the instanton 
solution for the tunneling through the potential barrier. 

Figure [7] is a contour plot of V(r) in units of V*, 
cf. Eq. d23j), for (3 > and u ± = w c /10, with r = 
(p, z) = r(sin$,cos$) (the angle ip is neglected due to 
the cylindrical symmetry of the problem). Darker re- 
gions correspond to a stronger repulsive potential. The 
repulsion due to the dipole-dipole and harmonic po- 
tentials is clearly distinguishable at |z|/7\t ~ and 
7, respectively. The lesser-dark regions located at 
(p±,±z±) = £±(sm$±, ± cosiJjJ correspond to the ex- 
istence of two saddle points, see circles in Fig. [7] Here, 
£±_ and cos^ arc £j_ = (^C^-o/mu^) 1 / 5 and cosi9±_ = 

\J\ — (r*/£j_) 3 /-s/5, respectively, while the barrier at the 
saddle point is V{p±,±z ± ) = C 3;0 /^i + C 6 . fi /£\. The 
figure shows that for distances r ^> £± > r+,rs the ef- 
fective interaction potential Eq. (|25p is purely repulsive. 
The existence of two saddle points at distances r ~ £± 
separating the long- from the short-distance regimes is a 



general feature of systems with (3 > and uj±/lu c < 1. 
Thus, £± defines the characteristic length-scale for at- 
taining purely repulsive 3D potentials in the presence of 
a static electric field. Actually, we show below that for 
collisional energies smaller thany(p_i_, \z±\) the dynamics 
of the particle can be reduced to a quasi two-dimensional 
(2D) one, by tracing over the fast particle motion in the 
z-direction. 

For strong trapping ui±_>ui c the two saddle points col- 
lapse into a single one located at z = 0, and p = £±_ ~ r+. 
In this limit the dynamics is purely 2D, with the particles 
strictly confined to the (z = 0)-plane. The long-distance 
regime is separated from the short-distance one by the 
potential barrier of height V(£ ± ,0) = K = S/3 4 / 54. The 
amount of harmonic confinement required to achieve this 
pure 2D regime increases rapidly with (3 as uj c oc /3 8 / 3 . 
While for a typical rotational constant B/h ~ 5 GHz 
and a weak DC field (3 = 1/10, oj c is of order of 
uj c /2tt ~ 10 kHz, for a (reasonable) electric field f3 = 1/3 
we find oj c /2it ~ 1 MHz. This value of cu c exceeds the 
tightest experimental optical traps w™ ax /27r ~ 150 kHz. 
Thus, in general the dynamics should be considered 
quasi 2D. 

When an ensemble of polar molecules is considered, 
inelastic collisions and three body recombination may 
lead the system to a potential instability, associated 
with the attractive character of the dipole-dipole inter- 
action [U IH, Hil ■ In our discussion, this instability is 
associated with the population of the short-distance re- 
gion r < £j_, which can be efficiently suppressed. In 
fact, for collisional energies smaller than the potential 
barrier V(p±, ±z_i_) the particles are mostly confined to 
the long-distance regime, where they scatter elastically. 
That is, when a cold ensemble of molecules is consid- 
ered the barrier provides for the stability of the system 
by "shielding" the short-distance attractive part of the 
two-body potential. In this limit, residual losses are due 
to the tunneling through the potential barrier. In the 
next section we estimate the tunneling rate T associated 
with this process, and we show that it can be efficiently 
suppressed for reasonable values of f3 and u>±. Thus, it 
is possible to realize stable 2D configurations of strongly 
interacting polar molecules interacting via dipole-dipole 
interactions 13011. 



4- Stability of long-range collisions 

In the following we calculate the rate T = Toe~ SE / h of 
particle tunneling through the barrier V(p±,±zj_) using 
a semi-classical/instanton approach [5l| . In particular, 
we focus on determining the quantity Se, the euclidian 
action of the semiclassical trajectory [5lJ, which is re- 
sponsible for the exponential suppression of the tunnel- 
ing. The constant Tq is related to the quantum fluctua- 
tions around the semiclassical trajectory, and its value is 
strongly system-dependent. For the crystalline phase of 
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Ref. [3(j (see also Fig. [3]), it is the collisional "attempt 
frequency" , proportional to the characteristic phonon fre- 
quency Tq ~ -J 'C^-q I ma 5 , with a the mean interparticle 
distance. 

The relative motion of the two particles in the effective 
potential V(r) of Eq. (f25|) is analogous to that of a single 
(fictitious) particle with reduced mass m/2, and dynam- 
ics determined by the Hamiltonian H = p 2 /m + V(r). 
The associated euclidian action, that is the action in 
imaginary time r, is given by 



Se [r(r) 



r+oo 


m 




/ dr 






J — oo 


7 





(26) 



where r(r) is the particle's trajectory. We remark that 
Eq. (|26|) corresponds to the action in real time, with an 
inverted potential —V(r). The classical trajectories are 
found by minimizing the action Eq. (|26p . yielding the 
following equation of motion 



m err 

2^ = +W «- 



(27) 



The "energy" E = p 2 /m — V(r) is conserved along each 
classical trajectory. The instanton solution is then the 
trajectory with the smallest action Se, which approaches 
r(r — > ±oo) = (p — > oo,0) asymptotically at time r — > 
=Foo, see dashed line in Fig. [7] The energy of the particle 
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FIG. 8: (color online) The euclidian action Se of Eq. (|28|) 
as a function of u±/uj c (solid line). For lj± < uj' c m 0.88 ui c 
(lj± > uj' c ) the "bounce" occurs for z(0) 7^ (within the 
plane z(0) = 0), see text. The point uj' c is signaled by a 
circle. The dashed line is the Ce-o- independent expression 
S E ~ 7.01S (w ± /tj c ) 1/5 (see text), with So = ^/7n\C^\/hrl. 
For iu± > lu' c the action is Se ~ 5.78So, which is ui±- 
independent, consistent with the "bounce" occurring in the 
(2 = 0)-plane. 



along this trajectory is zero. The action Se reads 

r(oc) 

'0 Jr(0) 

where r(0) is the "bouncing point" reached at r = [5lj . 



poo />r(oc) 

S E = 2 dr2V[r(T)} =2 dsy/mV(r), (28) 

JO Jr(0) 



We solve Eq. (|27|) numerically for the classical tra- 
jectories with zero energy, for a generic value of (3 and 
The obtained action Se is plotted in Fig. [7] as a 
function of iv±/u c , in units of So = ^/m|C 6;0 |/?* 2 = 
(2Sd 4 m 3 /3 8 /3 7 ) 1/6 . We notice that the action shows dif- 
ferent behaviors for uj±_ -C uj c and uj±_ 3> uj c . In particu- 
lar, for lo_\_ <C w c the action increases with increasing lv± , 
while for uj±_ lu c it is w^-independent. The transition 
between the two different regimes mirrors the change in 
the nature of the underlying potential V(r) as a func- 
tion of oj±/lu c described following Eq. (|25|) . as explained 
below. 

We find numerically that for u> < cu' c f=a 0.88 u c the 
"bouncing point" r(r = 0) = [p(0), z(Q)] of the instanton 
solution occurs for z(0) ^ 0, see dashed line in Fig. [7] 
This is consistent with the existence of two saddle points 
located at V(p±,±z±), with z± > 0. Since the saddle 
points appear approximately at a length r ~ £±_ 3> r+ , it 
is expected that in this regime the action is independent 
of the short distance behavior of the potential, that is 
of the C6;o-coefficient of the Van der Waals attraction. 
Accordingly, Fig. [7] shows that S*e is well approximated 
by S E « 7.01S (uj±/uj c ) 1 / 5 = 5.86(C 3 2 . m 3 w ± /8) 1 / 5 = 
lA3h(£±/a±) 2 (dotted line), which only depends on the 
C3 ; o-coefficient of the dipole-dipole interaction and the 
confinement along z, via a±_ = (h/muj^) 1 ^ 2 . 

For lu > lu' c we find numerically that the "bounce" 
takes place in the plane z = 0. This is consistent with 
the existence of a single saddle point located at V(£j_,0) 
for uj± > uj c , as discussed in the previous section. The 
"bouncing point" is at pir — 0) = l±_ — r^/2 1 / 3 and the 
action is S E = 5 2 5 / 3 v^r(7/6)/r(5/3) w 5.78 S - The 
latter is independent of loj_, which is again consistent 
with the collisional dynamics being purely 2D. 

The reason why the transition between the two behav- 
iors of the instanton solution happens at a value of uo±_ 
which is slightly different from uj c is that the instanton 
solution accounts for the kinetic energy of the fictitious 
particle. Thus, the particle is not always forced to follow 
the gradient of the potential, see dash-dotted (red) lines 
in Fig. [7] This results in the "bounce" occurring in the 
plane z = even for values of ujj_ slightly smaller than w c . 

From the discussion above it follows that in the limit 
of strong interactions and tight transverse confinement 
r rapidly tends to zero. We illustrate this for the ex- 
ample of SrO, which has a permanent dipole- moment 
of d as 8.9 Debye and mass m = 104 amu. Then, for 
a tight transverse optical lattice with harmonise oscil- 
lator frequency uj_ — 2ir x 150kHz and for a DC-field 
13 = dE DC /B '= 1/3 we have (C| ;0 m 3 w ± /8n 5 ) 1 / 5 as 3.39 
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FIG. 9: (color online) The ground-state effective 2D BO- 
potentials V^p(p) of Eq. ([29} as a function of the molecule 
separation p in the [z = 0)-plane for various strengths of 
the DC electric field (3 = 0,0.1,0.15,0.2,0.3,0.4. The quan- 
tity a± is the harmonic oscillator length in the z-direction. 
The molecular parameters are chosen as (d 4 m 3 B /ft 6 ) 1 / 2 = 
1.26 x 10 6 and the frequency of the harmonic potential in 
the z-direction is uj±_ = 15B/10 6 h. This corresponds to the 
case of SrO with a mass m = 104 amu, a rotational constant 
B w h 10 GHz and a permanent dipole moment d w 8.9 Debye 
in a tight confining potential with lu± — 2n x 150kHz, where 
a± ~ 25nm. 



and obtain r/r « e - 5 - 86x3 - 39 w 2 x 1CT 9 . Even 
for DC field as weak as j3 = dE-oc/B = 1/6 we still 
obtain a suppression by five order of magnitudes, as 

r/r « e - 5 - 86x1 - 94 « io- 5 . 



5. Effective 2-D interaction 

In the limit of strong interactions and tight optical con- 
finement, it is possible to derive effective two-dimensional 
potentials by integrating out the fast particle motion in 
the transverse direction z. 

For r > £j_ 3> ax, the two-particle eigenfunctions 
in the z-direction approximately factorize into prod- 
ucts of single-particle harmonic oscillator wave-functions 
ipk 1 (zi)i^k 2 { z 2)- In first order perturbation theory in 
V^g 3 '/fou)j_, the effective 2D interaction potential VJP 
reads 



dZidZ2M*l) 2 M*2) 2 Vrt(T) 

L- J dze- z2 /^V™{v). (29) 



Expression Eq. ([2"9"| is valid for large separations r > 
£j_ 3> ax where the potential is (much) smaller than 



the harmonic oscillator spacing, i.e. |VjP(p)| <C hui±. 
When this condition breaks down, more harmonic oscilla- 
tor states should be considered in addition to the ground 
states ipo(zi)tpo(z2) in deriving VJP from Eq. (|2"Tj) . In 
any case, for large separations p S> l± the 2D potential 
reduces to 



3;0 



<>:() 



Finally in the adiabatic approximation we obtain the 
effective 2D Hamiltonian H^g 



H, 



2D 
eff 



f-' 2m 



^: + K 2 ff D (p), 



(30) 



where = (p X ;j,Py,j) is the (two-dimensional) mo- 
mentum in the plane z = of molecule j — 1,2 and 
P = (%2 —Xi,y2 — yi) is the (two-dimensional) separation 
of the molecules in the plane z = 0. The derivation of 
H^g is the central result of Sec. IIII Al 



B. Effective interactions in the presence of an AC 
microwave field 



In this section we consider the interactions of two 
polar molecules in the presence of an AC microwave 
field of frequency u> and polarization q, with re- 
spect to the direction of transverse trapping e 2 , i.e. 
Eac(£) = EAce~ lu>t e q + c.c. The spatial dependence 
of Eac(£) is neglected, in accordance to the discus- 
sion of Eq. (J7J. The field is blue-detuned from the 
(Jj=0-1) transition of the single-particle rotor spec- 
trum by A = uj — 2B/h > 0, with Rabi- frequency 
n = E A c j (l,q\d qd \0,0) j /H = dE AC /V3h. 

The effects of the AC field on the two-particle scat- 
tering can be summarized as: (a) Inducing oscillating 
dipolc-moments in each molecule, which determine long- 
range dipolc-dipole interactions whose sign and angular 
dependence are given by the polarization g and the orien- 
tation in space, e r ; (b) Inducing a coupling of the ground 
and excited state manifolds of the two-particle spectrum 
at a resonant (Condon) point rc = (d 2 /3/iA) 1 / 3 , where 
the dipole-dipolc interaction becomes comparable to the 
detuning A. This coupling is responsible for an avoided 
crossing, whose properties depend crucially on the po- 
larization q. We show below that the character of the 
(3D) ground-state effective interaction potential is very 
different at distances larger and smaller than rc. 

The basic features of the scattering in the presence of 
the AC field are depicted in Fig. [TDJ In the figure, the 
solid (dashed) lines arc the bare (Eac = 0)-cigcnvalucs 
Ej-Y-a(r) of Eq. ([2D]) for a = + (a = -), plotted as 
a function of r. The color conventions arc the same as 
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in Fig. [5l The microwave field, which is detuned from 
the single-particle rotational spacing 2B by an amount 
TiA > 0, is represented by a black arrow. Analogous to 
Fig. [5j the excited-state manifold is split by the dipolc- 
dipole interaction. This splitting has the effect to ren- 
der the detuning position-dependent, so that eventually 
the combined energy of the bare ground-state plus a mi- 
crowave photon becomes degenerate with the energies of 
some bare excited states. The resonant points are de- 
noted as tq and r' c for the resonance with two symmetric 
(|$i ; i ±;+ (r)) = |$ l!+x . + (r)} ± |$i ; _i !+ (r)» and an anti- 
symmetric state (|$i ; o ; -(r))), respectively. The symmet- 
ric bare ground-state is coupled by the AC field to the 
symmetric bare excited-state |$i ; i_ ; +(r)) only, while the 
state |3>i ; i +; +(r)) is dark. As it is explained below, this 
coupling induces a splitting of the field-dressed energy 
levels at rc- Due to this coupling, the 3D effective dressed 
adiabatic ground-state interaction potential inherits the 
character of the bare ground and excited potentials for 
r > rc and r -C re respectively (thick solid line in the 
figure). Since the symmetric excited-state potential is re- 
pulsive, during the collision the dynamics of the particle 
is confined to the region r > rc, that is, the AC cou- 
pling can determine an effective "shielding" of the inner 
part of the molecular interaction potentials (the molec- 
ular core of Fig. [5]). This shielding is three-dimensional 
and it is analogous to the optical shielding of Napolitano, 
Weiner and Julienne developed in the context of ultra- 
cold atomic collisions [39l |40| . In particular, we show 
below that the shielding efficiency depends strongly on 
the chosen polarization q of the AC field (see Fig. [TT] and 
text), a characteristic which was found both in theory 
and in experiments with cold atoms [3^ . l4lj . 

Analogous to the optical shielding case, one expects 
that diabatic couplings among symmetric states provide 
for a loss mechanism in the 3D ground-state collision, 
for any finite collisional energies. In particular diabatic 
couplings, and therefore losses, are expected to be par- 
ticularly relevant in the region r ~ rc and r < rc, 
where the ground-state energy shows an avoided crossing 
with another potential, and the ground-state energy be- 
comes doubly degenerate, respectively (see Figs. [TU] and 
1111) . When a harmonic confinement in the z-dircction is 
considered, other loss channels may arise due to resid- 
ual non-compensated tensor-shifts Eq. (]12[1 . coupling 
the ground-state to the anti-symmetric state |$i ; o ; _(r)), 
whose energy Ei-o--(t) crosses the ground-state potential 
atr c = (2d 2 /3fiA) 1/3 f see Fig. [TP]). When more particles 
are considered, three-body interactions are expected to 
generate similar couplings to the anti-symmetric state. 
Three-body interactions are of concern since, as noted 
in Sec. Q] we are interested in designing effective two- 
dimensional interaction potentials for pairs of molecules, 
which can lead to the realization of interesting phases for 
an ensemble of polar molecules in the strongly interacting 
regime (see Fig. Q]and Ref. [3fj|). 

Because of all these loss mechanisms, two-dimensional 
shielding is not expected to be very efficient in the case 
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FIG. 10: (color online) Schematic representation of the ef- 
fects of an AC microwave field on the interaction of two 
molecules. The solid and dashed lines are the bare (Eac — 0)- 
potentials E n (r) = Ej-Y;cr(r) of Sec. MI A II for the sym- 
metric (a — +) and antisymmetric (a = — ) states, respec- 
tively. An AC field of frequency u> = 2B + A is blue de- 
tuned by A = 3B/10 6 h from the single-particle rotational 
spacing 2B, with Rabi-frequency fl. The dipole-dipole inter- 
action splits the excited-state manifold, making the detun- 
ing position-dependent. Eventually, the combined energy of 
the bare ground-state potential Eo-o-,+ (r) and of an AC pho- 
ton (vertical arrow) becomes degenerate with the energies of 
the bare symmetric -Ei ; ±i ; + (r) and antisymmetric Ei- t o--(r) 
potentials. The corresponding resonant points are denoted 
as r c = {d 2 /3HA) 1/:i (circles) and r' c = (2d 2 /3ftA) 1/3 , re- 
spectively. The resulting dressed ground-state potential is 
sketched by a thick solid line. For molecular parameters of 
SrO (B = h 10 GHz and d w 8.9 Debye) the detuning cor- 
responds to A/2tt = 30 kHz, and the lengths rs and rc are 
given by tb ~ llnm and rc ~ 0.5 fim, respectively. 



of interactions in an AC field. However, in Sec. IIII CI we 
show that most of these losses can be avoided, and an 
efficient 2D shielding recovered, by utilizing a properly 
chosen combination of static and microwave fields, and 
a tight harmonic confinement in the ^-direction. 

In the remainder of this section we further detail the 
interaction processes. This analysis is instrumental to the 
discussion of the collisions of two particles in the presence 
of both static and microwave fields, which is addressed 
in Sec. ImCl 

In Sec. IIII B II we derive the dressed adiabatic poten- 
tials for the interaction of two particles in an AC field. 
There, we show that the shielding is strongly dependent 
on the chosen polarization of the AC field. In fact, for lin- 
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ear polarization (q = 0) the width of the avoided crossing 
at the Condon point rc between the ground-state poten- 
tial E 0l0l+ (r) and the potential i?i ;+ i ;+ (r) is dependent 
on the value of the polar angle i9, and it vanishes for 
t> = and -d = it. This vanishing of the width of the 
avoided crossing entails the existence of "holes" in the 
three-dimensional shielding, which allow for reaching the 
molecular-core region [see Fig. [TTJal) and Fig. [TTT a2)]. 
On the other hand, for circular polarization (q = 1) the 
potential is repulsive in three-dimensions [see Fig.nTUbl) 
and Fig.[Tl](b2)]. Diabatic losses are most likely to occur 
at the Condon point rc , and for r < rc due to couplings 
to the dark state |$i ; i +;+ (r)), which becomes degenerate 
with the ground-state. In Sec. IIII B~2l and Sec. IIII Bifl thc 
interaction is further analyzed by deriving a perturbative 
expansion for the ground-state potential valid to second 
order in fi/A, and by analyzing a reduced model Hamil- 
tonian valid in the vicinity of rc, respectively. There, it is 
argued that couplings to the antisymmetric manifold due 
to three-body interactions (and to the possible existence 
of residual tensor shifts for a harmonic confinement) re- 
duce the efficiency of the shielding. 



1. Adiabatic potentials 

The total Hamiltonian for the collision of two particles 
in the presence of an AC field is 



(al) 



m = E 



i 



2m 



2 2 

mujj_Zj 



H- mt (r,t), (31) 



with 



fl int (r, t)=J2 [B3 2 j - (EAce-^dgj + h.c.)] + V dd (r) 



(32) 



Analogous to the discussion of Sec. IIII Al fl int (r,t) en- 
tirely determines the non-trivial dynamics of the sys- 
tem, since the harmonic motion of the center of mass 
is decoupled from the relative motion. The permuta- 
tion symmetry a — ± is conserved during the colli- 
sion, since fli n t(r, t) is invariant under the exchange of 
the position of the two molecules (j = 1) <-> (j = 2), 
i.e. r — > — r. Thus, H int (r, t) can be conveniently rewrit- 
ten as flint (r, t) = J2a=± p <? H int( r i t ) p <?i where P + and 
P_ denote the projector onto the symmetric and anti- 
symmetric manifolds, respectively. 

We obtain the solution of the time-dependent problem 
is obtained analogous to Sec. ?? by diagonalizing 
the Hamiltonian fli, lt (r, t) in a Floquet picture and 
proceeds as follows: First we diagonalize the Hamil- 
tonian in the absence of the AC field, Eac = 0, as 
^int(r) = Enl*n(r)>-Bn(r)(*„(r)| with n = (J;Y;a), 
which is the same as Eq. (|20[) and in particular is 
time-independent (see also Table HI|) . Then, we consider 
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FIG. 11: (color online) The dressed adiabatic potentials 
E n (r, &) of Eq. (|3T[) for two molecules interacting in an AC 
field. The setup is the same as in Fig. 1101 The field po- 
larization is linear (q — 0) in panels (al) and (a2), while it 
is circular (q = 1) in panels (bl) and (b2). The solid (red) 
and dashed (blue) lines correspond to the potentials for the 
symmetric (a — +) and antisymmetric (a = — ) states, re- 
spectively. The thick continuous (black) line is the adiabatic 
dressed ground-state potential _Ei ; o ; +(r). Panels (al) and (bl) 
show the potentials as a function of the separation r for in- 
teractions in the (z = 0)-plane (i9 = 7r/2). The position of 
the resonant Condon-point rc is indicated by a vertical line. 
Panels (a2) and (b2) show the angular dependence of the po- 
tentials at r = rc- Note that for q — (panel (a2)) Eo-o;+(r) 
becomes degenerate with £ , i ; i ±1 +(rc, #) at $ = 0,7T, while it 
is non-degenerate at all angles for q = 1 (panel (b2)), sug- 
gesting better shielding. The potential Ei-o-,- (rc, $) has an 
energy larger than Ei-o-,+ (rc, #) for all angles $, indicating a 
level crossing ar r > rc (see text). 



the effect of the AC field, EacW, vm a transformation 
to the Floquet picture, which is obtained by expanding 
the time-dependent wave-function in a Fourier series in 
the AC frequency to. After applying a rotating wave 
approximation, i.e. keeping only the energy conserving 
terms, we obtain the time-independent Hamiltonian 
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H(r), which describes the driven system. The Hamilto- 
nian preserves the permutation symmetry, a = ±, i.e. 

H mt (r) = J2 a Pa-Hint ( r )P<?- Analogous to the zero-field 
case of Sec. IIII A[ we restrict the basis set to the 16 
states belonging to the three lowest-energy manifolds. 
This is obtained by choosing a detuning much smaller 
than the rotational spacing A <C B, and working in 
the regime of weak saturation f2 <JC A. In fact, in this 
limit the anharmonicity of the single-particle rotational 



spectrum ensures that the population of high-energy 
rotational states is negligible. Finally, we solve for 
^int (r) by diagonalizing the Hamiltonian in the symmet- 
ric and antisymmetric subspaces separately, e.g. (r). 

The Hamiltonian (r) for the symmetric sub- 

space expressed in the basis |$j.y.+ (r)) with (J;Y) = 
{(0; 0), (1; Y)\ Y =-i,o,i, (2; lOI^-a.'-i.o+.o-.La} reads 
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(33) 



The Hamiltonian H^ nt \r) for the antisymmetric subspacc on the basis |$j ; y : _(r)) with (J',Y) 



{(1;*% 



(2;^)| 



F=-1,0,1 



}) reads 



■aH 








"0 















n*_ 


















n*_ 


— o* 









^2,-1 








-n+ 










A (_) 

^2,0 










O 








A (_) 



(34) 



In Eqs. (f33|) and (|34|) off-resonant couplings of order 
0(£lcP /r 3 B) between the various J-manifolds have been 
neglected. The detunings Ajy and couplings f2y in 
Eq. (|33|) and Eq. (p4]l depend on the separation r of the 
two molecules as 

= A$(r) = J A - Ej. Y;a (r)/h (35) 

n Y = si Y (d,<p) = siDl tY (<p,0,oy. (36) 

Here D\ >Y {<p,4,ti) = (1, g| exp(-»y>J r ) exp(-i#J y )\l, Y) 
are matrix-elements of rotation-operator, which ro- 
tates the lab-frame onto the frame where the colli- 
sion axis is fixed along eo. The coefficients c± are 
c± = [(l±l/V3)/2] 1 /2. 

As said above, a set of dressed BO-potentials E n {r) 
and of adiabatic eigenstates |$„(r)) is obtained by diag- 



onalizing the Hamiltonian H ln t as 

H iDt (r) = |$„(r))£ n (r)(# n (r)|, (37) 

n 

with n = (J;Y;a). The tilde refers to the implicit de- 
pendence of the dressed potentials and eigenstates on the 
Rabi-frcquency Q, the polarization q and the detuning A 
of the external AC-field. As mentioned above, we focus 
on blue detunings A = u> — 2B /K > 0, since we are in- 
terested in repulsive potentials which can "shield" the 
short-range molecular-core interaction. 

2. Asymptotic expansion: r>rc 

An insight into the nature of the dressed ground state 
potential can be obtained by deriving an expression for 
Eq(t) = i?o : o:+( r ) pcrturbativcly in the small parameter 
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CI/ A. The pcrturbativc expansion is valid at separations 
r tq = (d 2 /HA) 1 / 3 , where the dipole-dipole interac- 
tion in the first excited manifold n = (1; Y; a) is smaller 
than the detuning of the AC field. Then, to second order 
in CI/ A the dressed ground-state potential reads 



£o(r) 



-SAS(r) 



+1 



h E 



2\Cl A ($,y)\ 2 
^A«(r)-AW(r) 
2ft|f2| 2 2ft|ft| 2 d 2 (2-3q 2 ) 



A 



A 2 



3r 3 



■(1 - 3 cos 2 i?), 
(38) 



where terms of order 0(l/r 6 ) have been neglected. 
The first term in Eq. (|38|) describes a quadratic single- 
molecule AC-Stark shift, which is positive for blue de- 
tunings. The second term is understood as follows: The 
AC-field induces in each molecule an oscillating dipolc- 
moment of magnitude (dj) ~ dCle q / A, and on average 
the oscillating dipoles give rise to an effective dipole- 
dipole interaction in the ground state, which is propor- 
tional to oc (dj) 2 /r 3 . Equation (fS5|) shows that the over- 
all sign of the induced interaction can be changed by 
varying the polarization q. 

The pcrturbativc expression for the ground-state po- 
tential breaks down at r ~ re = (d 2 /3HA) 1 / 3 , where 
two of the bare (J = l)-excited potentials (-Ei ; ±i ; +(r)) 
become degenerate with the energy of the ground-state 
plus a photon of frequency w. 

While the validity of perturbation theory ceases at 
r ~ rc, further insight into the solution of the adiabatic 
scattering problem can be obtained by direct inspection 
of a specific example. Figure QT] shows the dressed 
BO-potcntials E n (r) of Eq. (J37J) for A = 3B/10 6 ft 
and Cl = A/4. The polarization is linear, q — 0, in 
panels (al, a2), while it is circular in panels (bl, b2), 
with q = +1. Panels (al) and (bl) show E n (r) as a 
function of the separation r for collisions in the plane 
i9 = arccos(z/r) = ir/2. Panels (a2) and (b2) depict 
the angular dependence of E n (r) at the Condon point 
r = rc for the two polarizations q = and q = 1, 
respectively. In all the panels, the solid and dashed 
lines denote symmetric and antisymmetric potentials, 
respectively, while the dressed ground-state potential 
Eq(t) = Et);0;+( r ) is represented by a thick solid line. 
Since we are interested in ground-state collisions, the 
figure suggests the two following observations: First, 
while the potentials Eq(y) is strongly repulsive for 
r < rc arid = tt/2, both for q = and 1 [panels (al) 
and (bl)], the angular dependence at r — rc is very 
different [panels (a2) and (b2)]. In particular, panel 
(a2) shows that for q = the repulsive potential is a 
maximum at i? = 7r/2, while it vanishes at •& = and 
$ = 7r. This vanishing of the repulsion allows for the 
molecules to approach the molecular-core region, and 
thus the polarization q = does not provide for an 
efficient three-dimensional shielding of the molccular- 
corc region. On the other hand, panel (b2) shows that 



the shielding may in principle work for q = 1, since 
the ground-state potential is repulsive for any angles. 
The second observation is that a level crossing of the 
ground-state potential with the antisymmetric potential 
i?i.o,-(r) appears at r' c = 2 1 / 3 rc for all polarizations 
[panels (al) and (bl)]. Couplings to this state can arise 
due to non-compensated residual tensor-shifts, when a 
harmonic confinement in the z-direction is considered, 
or due to three-body interactions, when an ensemble of 
polar molecules is considered. These couplings will in- 
duce losses in the ground-state interaction. In Sec. IIII CI 
we show that the position of this level crossing can be 
shifted to distances r < rc, and the associated losses 
can be avoided, by superimposing a weak static electric 
field to the AC field. In this way, an efficient (2D) 
shielding of the molecular-core region can be recovered. 



3. Resonant Condon-point: r ~ rc 

In the remainder of this section we analyze fur- 
ther the scattering process at the resonance point rc- 
We restrict the discussion to the three relevant states 
{|$0;0;+(r)),|$i ;+ i; + (r)),|$i ; _i ;+ (r))}, since all other 
symmetric states of the (J = l)-manifold are detuned 
by Aj CT y(r) > A 3> Cl. In this subspace the Hamiltonian 
Eq. 133]) reads 



-Hint(r) 



Ag(r) fi_ 



(0 t <p) A<+>(r) 



Cl 

fi+(0,¥>) 











(39) 



For q — we have Cl±(d, <p) = ^fisini? / 'y/2 and 
the ground-state couples only to the bright super- 
position state, |$i ; _i ; +(r)) — |$i ; +i : _(r)), with cou- 
pling v^20sin#. The orthogonal state, |<I>i ; _i ;+ (r)) + 
|$i;+i;— (r)), is dark with respect to the AC coupling. 
The dressed ground-state is then a (position-dependent) 
superposition of bare ground and excited states and the 
corresponding dressed potential is 



E Q , , + {v)/n= -A+(r) + V /A_(r) 2 + 2|r!| 2 sin 2 t ?, (40) 

with A±(r) = [A^ } (r) ± A ( ^(r)]/2 « A - d 2 /3hr 3 . 
We notice that for q = and $ = 7r/2 the potential is 
repulsive with a ~ 1/r 3 radial dependence, due to the 
avoided crossing at r = rc, see Fig fTTT al). For $ ^ tt/2 
the splitting of the avoided crossing decreases as sin $ and 
vanishes at r = rc for # = and $ = ir, sec Fig. [TTT a2). 
Thus, close to the point r = rceo the molecules can 
penetrate the 3D "shield" provided by the AC field and 
approach the short-range molecular-core region, rCrc- 
This behavior resembles the one encountered in Sec. IIII Al 
for the collision of two dipoles polarized by a DC field, 
when the intcrmolccular axis is parallel to the direction of 
the DC field, see Fig.[5](e). However, this unstable region 
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now appears at distances r ~ rc ~ (d 2 /3HA) 1 / 3 , which 
are larger than the short-distance length = (d 2 /B) 1 / 3 
by a factor ~ (B /HA) 1 / 3 . For a detuning on the order of 
tens of kHz and a rotational spacing of tens of GHz, rc 
is two orders of magnitudes larger than tb- 

For a circularly polarized field |g| = 1, we 
have fl±(i9,ip) = e lqv (l qp <?cos?9)/2 and hence the 
ground-state couples to the bright superposition state, 
cos 2 (^/2)|$i ;?;+ (r)) + sin 2 (^/2)|$i ;9;+ (r)), with an am- 
plitude \/2Q which is now independent of the angle 
The orthogonal superposition, cos 2 (i9/2)|$i ;9;+ (r)) — 
sin 2 (^/2)|$i ;g; +(r)), is dark with respect to the AC- 
coupling. The dressed ground-state potential is 

£ ;0;+(r)/ft = -A+(r) + v/A_(r) 2 + 2|r!| 2 , (41) 

with A± (r) defined as in Eq. ([4*0]) . The behavior of the 
ground-state potential in the (i? = 0)-plane is analo- 
gous to the linearly polarized case, see Fig. [TTT al.bl). 
However, in contrast to the q = case, now the width 
of the avoided crossing remains finite at all angles [see 
Fig- rn7 b2)]. and the AC shielding of the molecular core 
is effective. 

However, as noted above, the pure AC shielding mech- 
anism has an intrinsic flaw that limits its utility, once 
an ensemble of polar molecules is considered. The anti- 
symmetric state |<i>i ; o,-(r)) is strongly repulsive with en- 
ergy — HA\ q (r) f» 2d 2 /3r 3 — HA and thus gives rises to 

a real crossing at r' c — (2d 2 /3HA) 1 / 3 = 2 1 / 3 rc (see dot- 
ted lines in Fig. Illj) . This crossing at distances larger 
than rc is expected to give rise to (strong) collisional 
losses when an ensemble of polar molecules is considered. 
In fact, three-body interactions can couple the ground- 
state to the antisymmetric |i'i j o,-(r))-state. In addition, 
analogous couplings can be provided by residual non- 
compensated tensor shifts, when a harmonic confinement 
in the z-direction is considered. In the next section we 
explain how some of these problems can be circumvented 
by introducing an additional static electric field. In that 
case, an efficient and collisionally stable 2D shielding of 
the inner part of the potential can be recovered. 

C. Effective interactions in the presence of both a 
DC and an AC fields 

In this section we consider the interactions of two po- 
lar molecules in the presence of both a weak DC field 
Edc = Enceo with f3 = dEuc/B <C 1 and of an AC 
microwave field Eac(*) = Ej^ce~ luJt e q + c.c, where the 
polarization q is defined with respect to the z-direction, 
lj is the frequency, A is the detuning from the single- 
particle resonance, and O is the Rabi frequency. 

As explained in Sec. Ill BH the DC field partially splits 
the three-fold degeneracy of the ( Jj = l)-manifold of each 
molecule by an amount ~ HS — id 2 E^ c /2QB (the mod- 
ulus of the projection \M\ is conserved). When the AC 
field is superimposed to the weak DC field, this splitting 



can yield significative advantages regarding the stabil- 
ity of the ground-state collision: (a) The ground-state 
can couple to a single non-degenerate excited state of 
the two-particle spectrum, thus avoiding diabatic losses 
due to the presence of symmetric (dark) states close to 
the ground-state for r < rc ~ (d 2 /3HA) 1 / 3 (sec discus- 
sion in Sec. IIII Bp . In fact, because of the splitting HS, 
the energies of other symmetric states become compara- 
ble to the dressed ground-state energy only at distances 
r ~ rg = (d 2 /hS) 1 ^ 3 <C rc, where the dipole-dipole in- 
teraction becomes of the order of the splitting HS. (b) 
The location r' c of the real crossing of Sec. IIII Til is also 
shifted to small distances r < rg -C rc, thus suppress- 
ing losses due to three-body-induced (or residual-tensor- 
shift-induced due to transverse confinement) couplings to 
the ground-state. 

Both of the outlined processes are shown in Fig. [T27a) 
and Fig. [T2Tb). which depict the bare (Eac = 0) energy 
levels of the two-particle eigenstates with J < 2 as a func- 
tion of the distance r, for j3 = 1/10. The polarization of 
the AC-field is q = and its frequency u> is blue-detuned 
from the (|0o,o) ~~ * |</>i,o))-transition of the single-particle 
spectrum by an amount A = u — (uJ + 25/3) > 0. In 
the figure, the continuous (red) and dashed (blue) 
curves are the bare BO-potentials for the symmetric and 
antisymmetric states, respectively. The presence of the 
AC-field is signaled by a black arrow at the resonant 
(Condon) point rc ~ (d 2 /3HA) 1 / 3 . Analogous to the 
case of Fig. HJb), Fig. [T27 a) shows that the (J = 1)- 
manifold is asymptotically split by the DC-Stark-shift 
HS. As said above, for r ^> rg this splitting suppresses 
the coupling among states of the (J = l)-manifold due 
to dipole-dipole interactions. Moreover, we note that 
since the characteristic length rg is such that r$ <C rc 
and r' c < rg (not shown in the figure), the presence 
of the splitting ensures that rc 3> r c , as opposed to 
the (-Edc = 0)-case of the previous section Sec. IIII Bl 
As a consequence, for r > rj > r' c we expect diabatic 
and three-body-induced losses to be largely suppressed. 
Then, we show below that a strong optical confinement 
in the ^-direction allows for the realization of stable 
2D collision setups, analogous to Sec. IIII Al However, 
at variance with the DC case of Sec. IIII Al utilizing 
a combination of DC and AC fields allows for much 
greater flexibility in designing intcrparticlc interactions. 
In particular, we here focus on the realization of a 
2D potential achievable with a single AC field, whose 
character is very different at distances larger and smaller 
than rc- 

In the remainder of this section, we discuss further the 
above-mentioned processes. In Sec. IIII CTI we derive the 
3D dressed adiabatic potentials for interactions in the 
presence of combined DC and AC fields. In Sects. HTl C 21 
and IIII C 31 we illustrate the main features of the two- 
particle interaction, by specializing to the case where the 
AC-field polarization is linear (q = 0), and by solving a 
model Hamiltonian comprising only a limited number of 
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FIG. 12: (color online) (a) Schematic representation of the effects of a DC and an AC microwave fields on the interaction 
of two molecules. The solid and dashed lines are the bare potentials E n (r) = Ej-M;a{r,'ff) of Sec. MI A 21 with # = 7r/2 for 
interactions in the presence of the DC field only, for the symmetric (a = +) and antisymmetric (<r = — ) states, respectively. 
The DC field induces a splitting hS of the first-excited manifold of the two-particle spectrum. A microwave-field of frequency 
uj — ZJ + 25/3 + A is blue detuned by A > from the single-particle rotational resonance. The dipole-dipole interaction 
further splits the excited-state manifold, making the detuning space-dependent. Eventually, the combined energy of the bare 
ground-state potential _Eo;0;+(r) and of an AC photon (black arrow) becomes degenerate with the energy of the bare symmetric 
_Bi ; o ; + (r, 7r/2). The resonant point rc = (d 2 /3ftA) 1//3 occurs at r ~ 46 rs- A second resonant Condon point occurs at (much) 
shorter distances r c i < rg — (d 2 /hS) 1 ^ 3 with an anti-symmetric potential (not shown), (b) Blow-up of the potentials of panel 



(a) with M = (see text in Sec. IIII C2[) . The 



ground-state potential is sketched by a thick solid line. 



states, whose energies are close to the one of the ground- 
state. An expression for the ground-state interaction po- 
tential is obtained which shows that at large distances 
r ^ rc the potential has a behavior ~ 1/r 3 similar to 
the one obtained for two molecules in a DC field. How- 
ever now the effective dipolar strength is given by the 
combination of both the DC and the AC fields, and it 
can be much weaker than for r < rc. Thus, 3D inter- 
action potentials can be engineered that have a marked 
"step-like" character, being strongly and weakly repul- 
sive at distances smaller and larger than rc , respectively. 
Analogous to the DC case of Sec. IIII A[ an effective 2D 
interaction potential shielding of the short-range region 
r « rj « rc is obtained by adding a harmonic confine- 
ment in the z-dircction and tracing over the fast particle 
motion along z, (see Sec. IIII C~4l and Sec. IIII C~5|) . 



1. Adiabatic potentials 

The total Hamiltonian including the couplings to DC 
and AC fields reads 



m = E 



Pj 1 2 2 

2m 2 - 1 ■> 



tfintM), (42) 



with 



flint (*,*) = E i Bj2 j - ^ocdo, 



- (E Ac e- wt d q . 3 +\i. C .)} +^ dd (r). 

(43) 

Similar to the discussion following Eq. (fTS)) . the non- 
trivial system dynamics is determined by the dynamics of 
the relative degrees of freedom, which decouple from the 
harmonic motion of the center of mass. In this section 
we set uj± = in Eq. (|42"]) . and thus diagonalizing the 
Hamiltonian for the relative coordinates in the adiabatic 
limit corresponds to diagonalizing Hi nt (r,t). The case 
lu±_ =/= is treated in the next sections. 
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The Hamiltonian H- m t(r, t) is invariant under the per- 
mutation of the two molecules, (J = 1) <-» (j = 2), 
and thus it can be conveniently rewritten as iJi nt (r,t) = 

22<t=± P&Hint ( r ->t)Pv Here P + and P_ are the projec- 
tors onto the manifold of symmetric and antisymmetric 
states, respectively. Since (several) external fields are 
present, parity is not conserved. 

In Sec. IIII A 21 we have already diagonalized 
7?int(r, t) in the absence of the AC field, that is 
#int(r) = £J$n(r))£„(r)($„(r)|. The respective 
adiabatic potentials E n (r) are shown in Fig. [51 together 
with the corresponding quantum numbers n — (J; M; a). 
We remark that J is not a good quantum-number since 
the electric field breaks the parity for each molecule; 
Thus J = Ji + J2 merely indicates the asymptotic 
manifold. The corresponding adiabatic potentials and 
cigcnstatcs for the ground-state n = (0; 0; +) [valid for 
r tb = (d 2 j 'B) 1 / 3 ] and the lowest excited states [valid 
for r ^> rg = (d 2 /M) 1 / 3 ] are given in Tab. IIIII Our goal 
in this section is to extend that treatment to account for 
the driving by the AC microwave field, which we assume 
to be near-resonant with the transition from the ground 
to the first-excited manifold, i.e. u> ~ T3 + (2/3 — q 2 )S for 
polarizations q = 0,±1, see Fig. 2] The average energy 
separation hZu is defined in Eq. ([6]). 



Since for /3 < 1 the single-particle rotor spectrum is 
strongly anharmonic and the AC-field is near-resonant 
with the (Jj = <-> l)-transition, we restrict our dis- 
cussion to the rotor states with Jj = 0, 1 for each 
molecule, that is we consider 16 two-particle states. 
Moreover, we focus on the region r r$, where the 
dipole-dipole interaction is (much) weaker than the DC- 
ficld-induccd splitting hS in the excited states. There- 
fore, (up to corrections of order ~ d 2 /hSr 3 and ~ 
d 2 /Br 3 ) the states |3>j ; Af ;CT (r)) are given by the states 

I $ j;M;<x = \®J;M„;<r{r -> oo,fi, tp)), which are re- 
ported in Tab. IIIII These states are independent of r, 
that is, they depend only on the orientation of the two 
molecules. In analogy to the treatment of Sec. IIII Bl 
we utilize the approximate states \&j-M;tr(r)) to diag- 
onalize the time- dependent Hamiltonian H ult (r,t) in a 
Floquet picture: We expand the time-dependent wave- 
function in a Fourier series in the AC frequency u>. Af- 
ter applying a rotating wave approximation, i.e. keeping 
only the energy conserving terms, we obtain the time- 
independent Hamiltonian -ff(r), which again preserves 
the permutation symmetry, a = ±. The Hamiltonian 

-^int ( r ) f° r the symmetric manifold (a = +) is ex- 
pressed on the basis {|$j. Mfi . (T=+ (r))} with (J;M M ) = 
{(0; 0), (1; 1 T )|; (1; 0); (2; 2 A i)| /1 l_, ,+; (2; 1 T ); (2; 0)} as 
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= i-^j. m (?",■$) = Joj — Ej-M^-aif, denote position-dependent detunings (for each rotational excitation), 
Cl± = fl±((p) arc orientation-dependent couplings, which are detailed below, and c± = c±($) 
depends on the polar angle The parameter £ is defined in the caption of Tab. IIIII 

The Hamiltonian H^~ t \r) for the antisymmetric manifold (a = — ) expressed on the basis {I'&j-j 
(J; M M ) = {(1; M; (1; 0); (2; 2); (2; 1 T )} reads 
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In Eq. (|4"4")l and Eq. (|4"S")) we neglected off-resonant 
(second-order) corrections ~ VIm^)^ /Mr 3 to the Rabi- 
frequency. The couplings f2_(y>), fio, ^-i-Cv 5 ) are given 
by n_( V ) = g/i^ACe^/V 5 . flo = (1 - q 2 )f E A C, 
ft+(<p) = qfiE A ce lqip /V2, respectively. Thus for lin- 
ear polarization (q = 0) one has f2o = and fi± = 0, 
while for circular polarization (\q\ = 1) = and 
Q+(<p) = ±Q-(ip) = fle ±lv /sqrt2 for q = ±1, respec- 
tively. 



2. Model Hamiltonian for A <C (5 

In the following we illustrate the main features of the 
scattering in the combined DC and AC fields, using the 
example of Fig. [T2J The Rabi-frequency f2 is chosen real 
and positive and in particular smaller than the detuning, 

<C A. Moreover, we choose A <C S since we want to 
address the potentials in regions where the two molecules 
in the first-excited two-particle manifold arc aligned by 
the DC field and not by the dipolc-dipole interaction, 
see Fig. [T2] The figure shows that for r 3> r$ the bare 
(Eac = 0)-states |$j ; M;«r(r)) with M ^ are largely 
detuned from resonance by an amount of order S 3> A. 
Thus, the discussion of Eq. (|44|) and Eq. ([45]) can be sim- 
plified by restricting the Hilbert space to the four states 
with M = only. For our basis-set, these are: The 
three symmetric states |$j ; o : +(r)) with J = 0,1,2 and 
the antisymmetric state |$i : o : -(r)). Then, Eq. (|44|) and 
Eq. |g5D reduce to 



in terms of complex cubic roots by 
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(46a) 
(46b) 



The position-dependence of the detunings = 
Aj'^r, i?) has the usual dipolar form Aj ff J(r, i?) ~ T/r 3 , 
with T = 1 - 3cos 2 ?9. Explicitly, we have A^(r,t?) = 
J A - [C 3;( . /; o;.)T/r 3 - C 6 . (/ . 0;CT )(tf)/r 6 ]/7*. The coefii- 
cients C^- n and C6 ;n (t?) are given in Tab. [TTT] for n = 
0, 5, 15 (symmetric states) and n = 6 (antisymmetric 
state), respectively. For the following discussion, it is 
important to notice that for a weak DC electric field 
(3 <C 1, the C3 ; „-coefficients C 3; (o ; o;+) = 9o ~ (d/3/3) 2 
and C3 ; (2 ; o;+) = <?I ~ (d/3/5) 2 are quite small, since they 
are suppressed by a factor ~ (3 2 . On the other hand, 
the coefficients C 3 .( 1: o : ±) = ffo.92 ± f± ~ ±gP/3 for states 
belonging to the first-excited manifold are as large as the 
bare dipolar coefficients, see Tab. U and Tab. IIIII 

By diagonalizing Eq. (|4"6"|) (a) we obtain the three 
dressed symmetric potentials Ej-q- + (t) (with J = 0, 1, 2) 
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.1:0 



+ (r) = ^ e ±2-//3 



2 V 27 4 



-ftA(r,i?), 

4fi 2 + E^M) 



1/3 



(47) 



where P = 40 2 + £, 7 [A£~J M) - A(r,tf)] 2 /2, 

Q = 2fi 2 [A $ (r, t?) - A(r, 0)] - (r, 0) - A(r, 0)] 

and A(r,t?) = ^jA^(r,i?)/3. The dressed potential 

for the antisymmetric state, i?i o _(r) = — ^A^. (r, #), 
is the same as the bare one. 



The dressed potentials Ej-o a (r,tf) are plotted in 
Fig. [TJ for A = 4ft = 3B/l6 6 'h, B = h 10 GHz, lin- 
ear polarization (q = 0) and /3 = 1/10. These parame- 
ters are the same as in Fig [T5J In particular, Fig. fT3Ta) 
shows Ej-fi-air,^ = it/2) as a function of the distance r, 
for molecules on the plane z = rcos^ = (i9 = 7r/2). 
Figure I13f b) is a three-dimensional representation of 
the potential-energy surfaces Ej ; o-,+ (r, #) = Ej.Q. + (p, z) 
for the three symmetric states J = 0,1,2, for finite 
transverse displacements z = rcos# (■# 7^ 7r/2), while 
Fig. WSi c) is the same as Fig. fT3Tb). with the addition of 
the potential Ei-o-—(p, z) for the antisymmetric state. 

In Fig. [T3] (a) the dressed ground-state potential 
Eo ; o ; +(r, tt/2) is the thick solid curve with largest en- 
ergy, which undergoes an avoided crossing with the po- 
tential £i;0;+(r,7r/2) at a distance rc ~ (d 2 /3ftA) 1 / 3 . 
The precise value of rc is derived below. The figure 
shows that the Condon point rc separates an inner re- 
gion r < rc where the ground-state potential is strongly 
repulsive E . a . + (r < r c ,7r/2) ~ C* 3 (r < r c )/r 3 , from 
an outer region r > rc where the potential is only 
weakly repulsive _E 0; o ; +(r > rc, 7r/2) ~ C^{r > rc)/r 3 
with (7 3 (r > rc) -C Cs(r < rc). This marked de- 
pendence of the potential strength on r is the realiza- 
tion of the "step-like" potential of Fig. [21 and it is 
due to the fact that the dressed ground-state inherits 
the character of the bare ground-state and of the bare 
state | ( f>i : o : +(r)) for r > rc and r < rc, respectively. 
Thus, we have C$(r > rc) ~ C , 3 ; (o ; 0;+) ~ (^/3) 2 and 
C 3 (r < rc) ~ C 3: (i : o : +) ~ d 2 /3. A harmonic confine- 
ment in the z-direction will be added in Sec. IIII C 41 to 
ensure the stability of the 2D interaction. 

Figure [TST b) is a three-dimensional representation of 
the dressed adiabatic potentials Ej-m-,+ {y) for the sym- 
metric states, plotted as a function of p = rsint? and 
z = r cos The blue and red regions correspond to re- 
pulsive and attractive potentials, respectively The thin 
gray lines are equipotential energy contours. For z = 
(p-axis) we recognize the case of Figure [T3] (a), where the 
symmetric ground-state potential has the largest energy. 
The position of the Condon point rc is indicated by an 
arrow. The avoided crossing between the ground-state 
potential i?o ; o;+( r ) an d the potential £i ; o ; +(r) observed 
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FIG. 13: (color online) Dressed adiabatic potentials Ej-M-,a(r) of Eq. (|47[1 f° r the interaction of two molecules polarized by a 
(weak) DC field Edc = i?Dce z with /3 = dEoc/B = 1/10 and dressed by an AC with linear polarization q — 0. The AC-field 
detuning and Rabi-frequency are A = 3B/10 6 ft and Q. = A/4, respectively. For a typical rotational spacing of B ~ h 10 GHz 
these numbers entail A/2n — 30kHz and Q./2n — 7.5kHz. Panel (a): Dressed adiabatic potentials E n (r) of Eq. (|47[) plotted 
as a function of r for z = (# = 7r/2). The potentials corresponding to symmetric (antisymmetric) states are given by solid 
(dashed) lines, and indicated by Ej-M-,a(r) for J = 0, 1, 2, M = 0, a = + ( J = 1, m = 0, a = — ). The dressed ground-state 
potential Eo-o;+{f) has the highest energy (thick solid line). The other potentials are asymptotically detuned by a multiple 
of A. The ground-state potential Eo-o;+{f) shows an avoided crossing with the symmetric potential _Ei ; o ; +(r) at r — rc- (b) 
Dressed adiabatic potentials Ej-M;a{r) = Ej-,M;tr(p, z) for the symmetric states (a = +) of Eq. (|47[1 plotted as a function of 
p — r sin ■& and z = r cos Blue regions correspond to a repulsive potential, red regions to an attractive potential. For z — 
(p-axis) we recognize the case of panel (a), where the symmetric ground-state potential has the largest energy. The position 
of the Condon point rc is indicated by a arrow. Accordingly, the avoided crossing with the potential Ei-o- + (r) observed at 
r — rc in panel (a) is now visible in transparency, below the upper layer. For \z\ > the potential Ei ; o-,+ becomes less and 
less repulsive. For \z\ > p/y/2 we have 1 — 3cos 2 ■& < and thus £a ; o ; +(r) is attractive and the Condon-point vanishes since 
the two states are off-resonant, (c) Dressed adiabatic potential as in panel (b), but also showing the antisymmetric states 
[a = — ). We see that the dressed potential E\.q.-(y) for the antisymmetric state with (1; 0; — ) is strongly attractive in the 
plane, i.e. for z = 0, which corresponds to the profile shown in panel (a). With increasingly separation \z\jr > the potential 
become less and less attractive. For \z\/p > l/v2 we have 3cos 2 i9 — 1 > and the potential becomes repulsive. Thereby a 
crossing between the asymmetric state and the ground-state appears at a (second) Condon- "point" (dashed line), as the 
two states are resonant, however due to the permutation symmetry do not couple. 



at r = tq for & = 7r/2 in panel (a) is now visible in trans- 
parency, below the upper layer. The figure shows that 
for \z\ > the potential Ei-o.+ (p, z) becomes less and 
less repulsive, and thus the Condon point rc = rc(i!})e r 
occurs at shorter distances (sec below, Eq. (|50[) ). For 
\z\ > /o/v2, we have 1 — 3 cos 2 ?? < and therefore 
£a ; o ; +(r) becomes attractive. Thus, the Condon-point 
vanishes since the combined energy of the bare ground- 
state plus a photon and the energy £/i ; o ; +(r) of the bare 
state |$i ; o ; +(r)) are not resonant [for the dependence 
of the bare potential i?i ; o ; +(r) on the angle see also 
Fig. [6^b) and Fig. Od)]. This vanishing of the avoided 
crossing for |z| > p/\2 corresponds to the formation of 
a "hole" in the 3D potential shielding the molccular-corc 
region, and it allows for the familiar attraction of dipolc- 
dipole interactions. The presence of this "hole" is rem- 
iniscent of the vanishing of the Condon-point at d = 
and $ = ir for the case of a linearly-polarized AC field in 
the absence of a DC field [see Fig. [TTIa2) hi Sec. EH]. 



However, here there are no dark states present at rc, 
due to the DC-field-induced splitting HS of the (J = 1)- 
manifold. This fact eliminates a significant non-adiabatic 
loss channel for ground-state interactions. 



From Fig. fT3T b) we see that a real crossing with 
the anti-symmetric state |$i ; o ; -(r)) takes place for 
\z\ > p/y2 at a second Condon point, which we denote 
as r' c . This is at variance with the case of Sec. IIII Bl 
[see Fig. [TT] (al)], where r' c was r' c > rc for all angles, 
thus opening loss channels due to three-body-induccd 
(or tensor-shift-induced, when a harmonic confinement 
along z is considered) couplings to the symmetric 
ground-state for any d. The exact position of the point 
r',-, is obtained in the next section. 
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3. Effective 3-D interaction potential 

In the following we are interested in the effective 3D 
potential V^P(r) for two molecules in their ground-state 
dressed by the external fields. In the absence of a trap 
(uj ± = 0) V c 3 ff D (r) reads 



V c f(r) ee E ;0;+(r, 0) - 2 [E Q .o + E' 0fi ] 



(48) 



where the terms in brackets are the Stark-shifts _Eo,o ~ 
-d 2 p 2 /6 and E' 0fi w M7 2 /A induced by the DC and AC 
electric field, respectively. 

At separations r>rc the effective potential resembles 
the dipolar potential for two dipolcs aligned along e z and 
in second order in the saturation amplitude fi/A is given 
by 



V c f(r) 



C3;(0;0; + )T" 20 2 (C 3; (1;0;+) 



C3;(0;0; + ))^ 



A 2 



(49) 



where T ee 1 — 3cos 2 i9, and terms of order 0(Vt A ) and 
0(l/r 6 ) have been neglected. The first term in Eq. |49|) 
is the dipole-dipole interaction for the two weakly polar- 
ized molecules induced by the DC field (see Sec. IIII A~2l) . 
while the second term is the familiar dipolc-dipolc inter- 
action induced by the coupling to the AC field. The pro- 
portionality factor (C3.( 1;0; +)— C3 ; ( ;0;+))2^ 2 /A 2 appears 
due to the competition of the oscillating dipole-momcnt 
(~ dfl/A) induced by the AC field with the permanent 
dipole-momcnt already present because of the DC field. 

The perturbative expression Eq. (|49p breaks down 
when the level spacing becomes comparable to the cou- 

f2. In particular, for 



pling, that is for |A^(r) - A 



,(+)| 

0;0 



A^q (rc) = Aq^,' (rc) an avoided crossing occurs between 

the potentials i?o ; o;+( r ) and i?i ; o ; +(r), which defines the 
resonant Condon point, rc ee rc(?9)e r , where e r is the 
intermolccular axis. The Condon distance rc(0) is pa- 
rameterized in terms of the polar angle r9 as 



rc(tf) 



C3;(l;0; + ) — C*3;(0;0;+) 



HA(1 - 3 COS 2 tf) 



1/3 



(50) 



For = 7r/2 (z/r = 0) the Condon point is attained 
at r c = r c (7r/2) = [(C 3;(1;0;+) - C 3;( o ; o ; +))AA] 1 /3 « 

[d 2 /3HA} 1 / 3 , see Fig. [T2] (a), which depends on the de- 
tuning A and the difference in the C 3; „ coefficients of 
the first excited state and the ground state (C 3; ( 1;0; +) — 
C 3; (o;0;+) ~ d 2 /3 for a weak DC-field (3 < 1). ' For 
d 7^ ir/2 (z/r ^ 0) the avoided crossing occurs at 
smaller separations rc($) < r c until it vanishes for 
cos 2 tf = (z/r) 2 = 1/3, see Fig.^b). 

The position of the point r' c ee r c (0)e r is deter- 
mined by the crossing between the dressed ground-state 
potential Eq-q. + (t) and the potential for the anti- 
symmetric state £i ; o ; _(r). As mentioned above [and 
shown in Fig. I13f b)] , this crossing occurs in the region 



\z\ > p/2 (cos 2 ?? > 1/3). The distance r c (z?) is given by 
r c (0) « [(C 3;( o ; o ;+ ) -C 3;(1;0; -))(3cos 2 ??- 1)/RA]V3. 

The discussion above suggests that an effective 2D in- 
teraction potential V 2 ^(p) with no losses due to cou- 
plings of the ground-state to other symmetric or antisym- 
metric states may be obtained for distances r > r^, by in- 
troducing a parabolic potential in the z-direction confin- 
ing the particles to the sector (z/r) 2 < 1/3. This "shield- 
ing" of the loss channels is analogous to the "shield- 
ing" of the attractive part of the potential and of the 
molecular-core region of the DC case for r > £± 3> r± ~ 
(d 2 /B(3 2 ) 1 / 3 (see Sec. IKTA"]) . However, now £ ± is ef- 
fectively replaced by rc ^ £±, and rc allows for much 
greater flexibility in tuning by external fields. 

In the next section we detail the requirements for ob- 
taining a stable effective interaction in 2D. In this way, 
it is possible to realize the 2D potential with "step-like" 
character, as shown in Fig. [5] 



4- Parabolic confinement 

The presence of a finite trapping potential of frequency 
luj_ in the z-direction provides for a position-dependent 
energy shift of Eq. (|49| . Thus, the new potential reads 



V(p,z)^V™(v) + \mulz 2 . 



(51) 



Analogous to the discussion of Sec. IIII A 3[ the combina- 
tion of the dipolc-dipolc interaction, which is repulsive for 
r > rc, and of the harmonic confinement yields a repul- 
sive potential which provides for a three-dimensional bar- 
rier separating the long-distance repulsive regime from 
the short-distance regime, where collisional losses can oc- 
cur. If the collisional energy is much smaller than this 
barrier, the relative motion of the particles is confined 
to the long-distance region, where the potential is purely 
repulsive. 

Figure (fT4")) is a contour plot of Eq. (f5Tj) for the same 
parameters as in Fig. [IjH i.e. A = 4f2 = 3B/10 6 h and 
(3 = 1/10. The frequency oj± for the harmonic confine- 
ment is uj± = A/5. In the figure, darker regions corre- 
spond to a stronger repulsive potential, and the white 
region for p,z < rc/2 corresponds to V(p, z) < 0. The 
repulsion due to the dipolc-dipolc and harmonic poten- 
tials is clearly distinguishable at z — and |z|/rc > 1, 
respectively. Two saddle points located at (p±, ±zj_) sep- 
arate the repulsive long-distance from the short-distance 
regions (circles in Fig. (jT4j) ) . The location of the saddle 
points approaches (px, ~ ( r 'Cj r c/2) with increas- 
ing confining potential uj± in the z-direction. In the fig- 
ure, the dotted line signaled by r' c marks the location of 
the crossing between the dressed ground-state potential 
Eo-a-+(r) and the potential Ei-o ; -(v) for the antisymmet- 
ric state. The figure shows that this crossing occurs in 
the short-distance region r < rc for all z, in agreement 
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in V^g 3 (r)/?ia;_L the 2D effective ground-state potential as 
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FIG. 14: (color online) Contour plot of the effective potential 
V(p, z) of Eq. (|51[1 for two polar molecules interacting in the 
presence of a weak DC field and an AC field. The field param- 
eters are the same as in Fig. [13] The frequency of the confining 
harmonic potential in the z-direction is lux = A/5. Darker 
regions represent stronger repulsive interactions. The white 
region for p < 1/2 indicates a potential V(p,z) < 0. The 
combination of the dipole-dipole interactions induced by the 
DC electric and AC (microwave) fields and of the quadratic 
confinement leads to realizing a 3D repulsive potential for 
r ^> to- In particular, two saddle-points located at (fix, ±z±) 
(circles) separate the repulsive long-distance regime r ^> rc 
from the short-distance regime r < rc where diabatic losses 
occur. The dotted line signaled by indicates the location 
of the crossing between the ground-state potential and the 
energy £i ; o ; -i(r) of the antisymmetric state, see text and 

Fig [H 



with previous discussions. Thus, for r > rc and colli- 
sional kinetic energies smaller than the potential barrier 
at the saddle point the ground-state interactions are sta- 
ble and purely repulsive, consistent with the discussion 
above. 



5. Effective 2-D Interaction 



Analogous to the discussion for the DC-field case in 
Sec. IIII A 51 in the limit of tight optical confinement it 
is possible to derive effective two-dimensional ground in- 
teraction potentials V^P(p) by integrating over the fast 
transverse degrees of freedom, z\ and z^. For r > rc > 
ax, the two-particle eigenfunctions in the z-direction ap- 
proximately factorize into products of single-particle har- 
monic oscillator wave- functions, and thus the integration 
is conveniently carried out in the harmonic oscillator ba- 
sis. In the adiabatic approximation wc find to first order 



V c f( P ) 



dz 1 dz 2 |Vo(^l)| 2 |V'o(^)| 2 V; 3 ff D (p, 2 2 - Zl ) 

1 



27raj_ 



dze-*-I^V™{p,z), 



(52) 



where ipk(zj) is the fc-th harmonic oscillator wave- 
function for the transverse confinement. In an analogous 
way, effective 2D potentials can be derived for all the 
dressed potentials Ej. M .,j(r), as (up to a constant shift) 



Ej-m 



i 



2ttqx 



dze- z2 l 2a ^E.j, M , a {p,z), (53) 



with So ; 0;+(p) = V cb(p)- In thc following we discuss 
the validity of thc adiabatic approximation in the case 
when both the DC and the AC fields are present. We 
focus only on thc four above-mentioned potentials, since 
the remaining states of the (J = l)-manifold are detuned 
by a large amount ~ 5 /lux ~ 10 3 . Thus, we neglect 
non-adiabatic couplings from the ground-state to the 
continuum corresponding to high-energy transverse ex- 
citations of the far-detuned states, since these couplings 
are expected to vanish at large inter-particle separations. 



At variance with the DC-field case of Sec. IIII A 51 sat- 
isfying the adiabatic approximation in the presence of 
both DC and AC fields is non as trivial. In fact, for a 
blue-detuned (A > 0) AC field the dressed ground-state 
potential -Eo ; o;+( r ) has the largest energy, see FigliUTa). 
Thus, it can happen that E . . + (r) becomes degenerate 
with the energy of one of the other states plus some mul- 
tiple k of the harmonic oscillator energy in the trans- 
verse direction Hlux- When these degeneracies happen, 
avoided and real crossings occur with the energies of the 
symmetric and the anti-symmetric states, so that satis- 
fying thc adiabatic requirement becomes in general much 
harder than in the DC-field case of Sec. IIII A 51 In fact, 
there the ground-state is the lowest-energy state and the 
lowest-energy excitations are the ones of the harmonic 
oscillator along z. In that case the adiabaticity condition 
is satisfied for VjP(p) <C haj±. On the other hand, it 
is still possible to derive an expression analogous to the 
latter even for the case when the AC field is present, if 
the trapping potential is large enough so that lux ^ 2A. 
In fact, then for large distances p > rc thc energy dif- 
ference between £?o;0;+(r) and the one of the first-excited 
state is approximately tkux — 2hA, where —2hA is the 
energy i? 2 ;0:+( r ) = — 2/iA (see Fig. [T3J) . In this case, the 
adiabatic approximation is still valid provided 



V c ff D (p) < fiwj. - 2fiA. 



(54) 



The perturbative expressions Eq. (|B"5j) for the dressed 
effective 2D potentials Ej B . a (p) are shown as thick 
dashed lines in Fig. [15] for the combination of a weak DC 
field with (3 = 1/10 and an AC field with linear polariza- 
tion q = and detuning A = 4£! = 3B/10 6 h. The panels 
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FIG. 15: (color online) Effective 2D potentials S3f M;(T;fc (p) with (J = 0, 1, 2; M = 0; a = ±1; k = 0, 1, 2, . . .) of Eq. JUJ for the 
AC- and DC-field setups of Fig. [HI The strength of the (weak) DC field is /3 = dE DC /B = 1/10. The AC field has polarization 
q = 0, is blue detuned by A = 3B/10 6 h and the saturation amplitude is fi/A = 1/4. Panels (a), (b) and (c) correspond to 
a harmonic oscillator frequency of the confining potential ui±/A = 1/2,2 and 5, respectively. The tick-dashed lines indicate 
the (single band) effective potentials Ejf . a (p) of Eq. (JSSJ). The solid lines are the non-perturbative potentials EjP . a . k (p) of 
Eq. (|55|l . where k = 0, 1, 2, . . . denotes the k th transversal excitations in the z-direction. The ground-state potential approaches 
twice the single-particle Stark-shifts, 2(Eo-o+E' . ) (see text), for p — > oo and it is indicated by arrow(s). The harmonic oscillator 
length a± for a typical mass m ~ lOOamu, a dipole moment d m 8.9Debye and a rotational constant B m h 10 GHz is indicated 
on the p-axis. The corresponding detuning is A w 2n x 30 kHz and the trapping frequency is u j_/2n w 15kHz, 60kHz, 150kHz 
in Panel (a), (b), (c), respectively. The gray region corresponds to r < r$- 



(a,b,c) of Fig. [T5] represent different transverse trapping 
frequencies given bya>^/A = l/2,2,5, respectively. The 
effective potential for the ground-state i?o ; 0;+ i s indicated 
at large separation p 3> tq , where it approaches the value 
i?o ; o;+(p — * °°) = 2i?o,o + 2i?o.oj corresponding to the 
DC and AC Stark-shift of the separated molecules, with 
2E' Q /hA ps +2{fl/A) 2 = 1/8. The thin solid and dotted 

lines in Fig. [T51 show the potentials Ej B M . a . k (p) for a = + 

and a = — , respectively. The potentials E 2 J ? M . cr . k (p) have 
been obtained by diagonalizing numerically 

^ci = ^ + ^n^ 2 -^ + ^i„t(r), (55) 
m 4 2 

with _ffj nt (r) given in Eq. ([46]) . Here, the index k = 
0,1,2,... labels the transverse excitations and at large 
separations p » re the corresponding potentials ap- 
proach ^j° /:cr . fe (p — > oo) ps Ej. M . a {r — > oo) + khw^. 

In Fig. [ToT a) we observe a series of avoided crossings 
involving the ground-state potential in the region p < tq. 
For p > rc the ground-state potential has the char- 
acteristic l/p 3 -dependence. Fig. IToTb) shows that for 
uj± = 2 A the ground-state potential V^g(p) is already 
well separated from the energy of the first-excited state 
with k > in a region p > rc . Finally, Fig. [T5T c) shows 
that for a tight trapping, lo±_ = 5A, the ground-state po- 



tential V^g(p) for p > rc/2 is well separated by HA 
from all the excited state with k > . The adiabatic ap- 
proximation is valid for V 2 s °(p) -C ftA, consistent with 
Eq. (H5). 

Remarkably, we find that since the spontaneous 
emission rates in the excited rotational levels of polar 
molecules, Tse, are negligible compared to achievable 
optical confinements ~ cu±_ ~ 27r x 150 kHz, the regime 
where A < oj±/2 is widely accessible. In fact one can 
achieve the limit of weak saturation and strong confine- 
ment while at the same time be sufficiently detuned to 
not suffer of spontaneous emission. That is, it is possible 
to fulfil all of the inequalities Tse <C Q < 2 A < iv±. 



IV. CONCLUSIONS 

In this work, we have shown how to engineer 2D inter- 
action potentials for optically-trapped polar molecules in 
their electronic and vibrational ground-state. In partic- 
ular, we have shown how to modify the shape as well as 
the strength of the inter-particle interaction potentials, 
by manipulating the rotational dynamics using external 
DC and AC microwave fields in combination with dipole- 
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dipolc interactions. In Sec. IIII Al we have shown that in 
the presence of a DC field and of a tight optical con- 
finement it is possible to realize effective 2D potentials, 
where particles interact via a purely repulsive ~ 1/r 3 . 
A potential barrier shields the attractive inner-region of 
the interaction potential, thus providing for the stabil- 
ity of the collisional setup. In Sec. IIII Bl we have ana- 
lyzed the interactions in the presence of an AC field. We 
have derived the 3D adiabatic potentials for the molecu- 
lar interactions, finding several degeneracies in the two- 
particle spectrum at distances of the order of the reso- 
nant Condon point re between the energy of the ground- 
state plus one photon and states of the first-excited man- 
ifold. The presence of these degeneracies opens (diabatic, 
three-body-induced, and, for the case of transverse con- 
finement, residual tensor-shift-induced) loss channels for 
the ground-state collisions, which make the case of in- 
teractions in a pure AC-field less appealing for realizing 
stable collisional setups in two-dimensions. In Sec. IIII CI 
we show that it is possible to realize stable 2D interaction 



setups with considerable flexibility in potential-designing 
by combining DC and AC fields, in the presence of strong 
transverse confinement. In fact, the DC field helps to 
greatly suppress the presence of loss channels at large 
distances, while the AC field allows for realizing poten- 
tials whose shape can vary markedly between the long 
and short distance regimes. 
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